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Abstract 



In this thesis we present contributions in the field of the apphcations of quantum field theo- 
ries techniques to condensed matter models. In chapter 3 we investigate on the non covariant 
fermionic determinant and its connection to Luttinger liquids. We address the problem of 
the regularization of the theory. In chapter 4 we treat spin hipping interactions in the non 
local Thirring model and we obtain an effective bosonic actions that describe separated spin 
and charge degrees of freedom. In chapter 4 we apply the self consistent harmonic approx- 
imation to previously derived bosonic action and we obtain potential depending equations 
for the spectrum gap. In chapter 5 we include spin-orbit couplings and compute correlations 
functions. We show that the spin orbit interactions modify the exponents and the phase 
diagram of the system and makes new susceptibilities diverge for low temperature. Finally 
in chapter 6 we summarize the main results and the conclusions. 
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Capitulo 1 
Introduccion 



La descripcion de sistemas dc fcrmioncs altamente corrclacionados es un problema central 
de la fisica de materia condensada. Durante las decadas pasadas, diversos experimentos 
Uevados a cabo en diferentes clases de materiales mostraron que la alta correlacion es un 
ingrediente fundamental a considerar para la comprension de sus propiedades fisicas. Entre 
otros podemos mencionar sistemas de efecto Hall, superconductores de alta Tc, y diversos 
metales, superconductores y aisladores organicos. Por otro lado, el tratamicnto tcorico de 
estos sistemas constituye una tarea formidable. Las ideas basicas acerca del comportamiento 
de Ids electrones en materiales son conocidas, al menos intuitivamente, desde hace muchos 
anos. En aparente contradiccion con lo expresado mas arriba, Sommerfeld [1] mostro que el 
comportamiento lineal del calor especffico de los metales a bajas tcmperaturas, al igual que 
el comportamiento asintotico a bajas temperaturas de la resistividad y de la conductividad 
optica podian ser entendidos suponiendo que los electrones en el metal se comportaban como 
un gas de fermiones no inter actuantes. Simultaneamente, Pauli [2] calculo la susceptibilidad 
paramagnetica de electrones litres y hallo que es independiente de la temperatura, en perfecto 
acuerdo con los experimentos. Al mismo tiempo, a partir de los trabajos de Bloch [3] y 
Wigner [4], se encontro que las energfas de interaccion de los electrones en el rango metalico 
de densidades era comparable a la energia cinetica. 

La resolucion de esta paradoja surgio con los trabajos de Landau [5,6] donde se intro- 
dujeron las ideas fundamentals que dominarian la vision de los sistemas interactuantes en 
materia condensada hasta nuestros dias. Landau postulo que los sistemas en interaccion 
evolucionan a partir de los sistemas libres al conectar la interaccion de manera adiabatica. 
Y que los estados de particulas en el sistema no interactuante se corresponden uno a uno 
con estados de cuasiparticulas o excitaciones elementales en el sistema en interaccion, es 
decir que poseen los mismos niimeros cuanticos. Existen sin embargo algunas restricciones: 
la mas importante impone que solo pueden considerarsc excitaciones en una escala dc en- 
ergia pequena comparada con la energia de Fermi. Esta restriccion sin embargo nada dice 
acerca de la intensidad de las interacciones que ocurren entre los electrones; estas pueden 
ser arbitrariamente fuertes. De alli que los sistemas fuertemente interactuantes que pueden 
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tratarse con esta teon'a se comportan cualitativamente como sistemas libres, con pardmetros 
renormalizados por las interacciones. Sin embargo, las propiedades de baja temperatura de 
muchos materiales que exhiben este tipo de comportamiento, poseen coeficientes que difieren 
hasta un factor de 10'^ respecto de los valores para electrones libres. Por otro lado los esta- 
dos en el sistema libre e interactuante deben tener la misma simetrfa, y ademas, al conectar 
la interaccion no deben formarse estados ligados. Estas restricciones impiden atacar con 
este formalismo problemas tales como el ferromagnetismo o la superconductividad, que se 
caracterizan justamente por esos efectos. 

La teon'a microscopica que respaldo esta teon'a fenomenologica pronto se desarrollo a 
partir de un Hamiltoniano fermionico con interacciones de dos cuerpos [7]. En general 
los objetos de interes son las funciones de Green G{k,u); si se las conoce para todos los 
valores de k y u;, pueden obtenerse en principio todas las propiedades termodinamicas del 
sistema. Su comportamiento a bajas energi'as y grandes longitudes de onda esta relacionado 
con el estado fundamental y los estados excitados mas bajos, y como el espectro de baja 
energi'a esta determinado cualitativamente por unos pocos parametros universales como la 
dimension, simetrfas y leyes de conservacion, el comportamiento infrarojo de las funciones de 
Green permite efectuar una clasificacion de los sistemas de muchos cuerpos en interaccion. 

En la mayon'a de las situaciones de interes es imposible calcular las funciones de Green 
en forma exacta, de modo que es necesario recurrir a metodos aproximados. El enfoque 
usual consiste en hacer un desarrollo perturbativo de (^(k, cj) en potencias de la interaccion. 
En los llamados Uquidos de Landau Fermi, este enfoque perturbativo es posible, y aunque 
para interacciones fuertes deben sumarse infinitos ordenes del desarrollo, las integrales gen- 
eradas en la expansion perturbativa estan libres de divergencias. Este desarrollo arroja como 
result ado aproximado 

G(k, uo + zO+) ^ ^ , (1.1) 

para la funcion de Green retardada con k en la vecindad de la superficie de Fermi. El mimero 
Zk es el llamado residue de la cuasiparticula, y la energfa es la energi'a de excitacion de 
una cuasiparti'cula. Como en general los Ifquidos de Landau Fermi son metales, ^k no debe 
tener gap. Esto significa que existe una superficie en el espacio k en la que ^k = 0, lo que 
define la superficie de Fermi, y Zkp representa la magnitud del salto de la distribucion de 
momentos en dicha superficie. La energi'a 7k puede identificarse con el amortiguamiento de 
la cuasi-parti'cula (o lo que es lo mismo, = l/7k con su tiempo de vida). Notese que 
en el piano complejo u, G(k,u + zO"*") posee un polo simple en u; = ,^k — ^7k con residuo 
Zk. La funcion de Green de los electrones no interactuantes, denotada Go(k,Lj), puede 
obtenerse como caso especial de la Ec. ()1.1|1 tomando .Zk = 1, 7k = e identificando ^k con 
la dispersion del sistema sin interacciones. En este caso el polo simple en = ^k — ^O"*" con 
residuo unidad representa la propagacion no amortiguada de una partfcula con energfa ^k- 
El correspondiente polo en la funcion de Green del liquido de Fermi en interaccion se asocia 
con el llamado polo de la cuasiparticula. El punto importante es que en la vecindad del 
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polo de la cuasiparticula, las funciones de Green del sistema interactuante poseen la misma 
estructura que las del sistema libre. Si definimos la autoenergia E(k, a;) como 

E(k, uj) = [Go(k, uj)]-' - [G(k, uj)]-' , (1.2) 

entonces las cantidades Z^, y 7k pueden calcularsc dc las derivadas de la autoenergia. 

En algunos casos, sin embargo, la aplicacion de esta maquinaria conduce a integrales 
divergentes en la expansion perturbativa de S(k, a;). El colapso del desarroUo perturbativo 
es un indicador de que la funcion de Green de la teon'a en interaccion no se relaciona mas 
en forma simple con la funcion de Green de la teon'a libre, por ejemplo por la existencia de 
polos multiples o de singularidades no algebraicas, lo que impide definir las cuasiparticulas. 
En este caso el sistema no puede ser un liquido de Fermi. 

En alios mas recientes, el interes en sistemas en una dimension espacial se vio incrementa- 
do debido a la realizacion experimental de materiales en los cuales el movimiento electronico 
correlacionado se encuentra efectivamente confinado a una dimension. Podemos citar como 
cjcmplos nanotubos de carbono [8], cstados de bordc cn sistemas de efecto Hall [9-11], con- 
ductores organicos [12], heteroestructuras semiconductoras (alambres cuanticos) [13], etc. 
Al intentar aplicar la teoria del liquido de Fermi a tales sistemas, se arriba a los problemas 
antes mencionados. El modelo mas sencillo considerado con el objeto de describir el estado 
metalico normal de estos sistemas es el modelo de Tomonaga-Luttinger [14-16]. Como vere- 
mos en el capftulo siguiente, este modelo es exactamente soluble, sus funciones de correlacion 
pueden ser calculadas, y todas sus propiedades se vuelven accesibles, por ejemplo el espectro 
de bajas energias resulta lineal y sin gap. Las excitaciones fundamentales no son mas las 
cuasiparticulas, sino fiuctuaciones bosonicas colectivas independientes de grados de libertad 
de carga y spin. Esta independencia entre ambos se denomina separacion spin-carga [17]. 

El colapso de la teoria de Landau del liquido de Fermi en este modelo pucdc cntcndcrsc al 
observar el comport amiento de su funcion dc Green, que presenta un decaimiento algebraico 
no universal, es decir, el exponente del decaimiento dcpcndc dc las intcraccioncs. Adcmas 
el calculo de la distribucion de momentos tambicn arroja un comportamicnto algebraico en 
las proximidades de la superficie de Fermi, y se vuelve continuo en dicha superficie, aiin a 
temperatura cero. 

Haldane [18] conjeturo que este conjunto de propiedades no es exclusivo del modelo de 
Tomonaga-Luttinger, sino que son propiedades genericas del estado metalico normal de sis- 
temas de electrones interactuantes en una dimension. Mas aiin, Uamo Liquidos de Luttinger a 
estos sistemas, y propuso que el modelo de Tomonaga-Luttinger es su Kmite a bajas energias, 
en el mismo sentido en el que el gas de Fermi es el modelo libre sobre el que se construye 
el liquido de Fermi. La inclusion de interacciones que sacan al sistema del punto fijo del 
liquido de Luttinger, como dispersion hacia atras o umklapp generan gaps en los espectros 
de carga o spin, y dependiendo de los valores de los acoplamientos pueden ocurrir transi- 
ciones de fase entre un estado sin gap de tipo Liquido de Luttinger y otros con gaps en los 
espectros de carga o spin, por ejemplo la Uamada transicion metal aislador de Mott [17]. La 
importancia de los liquidos de Luttinger recobro I'mpetu hace pocos aiios desde la propuesta 
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de Anderson de que poseen propiedades que son semej antes a la de los super conductores de 
alta temperatura crftica [19,20]. 

La correspondencia entre teoria y experimento para estos sistemas unidimensionales se 
ha visto puesta a prueba en diversos materiales. Las leyes de potencia de funciones de 
correlacion se han verificado en muchos de ellos, por ejemplo en las Ref. 21-23 mediante 
propiedades de transporte. Sin embargo la separacion spin-carga ha sido mas elusiva, y las 
pruebas mas convincentes son muy recientes [24]. 

El marco teorico en el cual se da tratamiento a estos sistemas se encuentra mtimamente 
vinculado a los formalismos utilizados en teoria cuantica de campos y en el tratamiento de 
las interacciones fundamentales de la naturaleza. La bosonizacion abeliana, o trasmutacion 
de fermiones en bosones, es la tecnica construida para la resolucion original del modelo de 
Tomonaga-Luttinger en materia condensada [16,25,26], y fue desarroUada en paralelo con 
la tecnica del mismo nombre en teoria de campos [27-29]. Esta se aplico a la resolucion del 
modelo de Thirring, y el modelo de Thirring masivo, que son versiones en el lenguaje de la 
teoria de campos del modelo de Tomonaga-Luttinger. Mas tarde, a partir del trabajo de 
Fujikawa [30] , la tecnica de bosonizacion en teoria de campos se extendio para ser aplicada en 
el contexto de las integrales funcionales [31,32], y un proceso equivalente ocurrio en materia 
condensada [33,34]. Mas especificamente, la bosonizacion funcional se aplico en la Ref. 35 
al estudio de una version del modelo de Thirring, en la que se modifico el acoplamiento entre 
las corrientes fermionicas para dar lugar a la posibilidad de una interaccion no local. En 
la Ref. 36 se aplico dicha formulacion al calculo de las funciones de Green del modelo con 
interacciones de largo alcance, mas especificamente de tipo columbiano. 

En esta tesis presentamos contribuciones originales en el campo de las aplicaciones de 
las teon'as cuanticas de campos a la formulacion de modelos de materia condensada en los 
que el spin electronico juega un rol crucial. En particular nos concentramos en el modelo de 
Tomonaga-Luttinger con spin y dos extensiones posibles: interacciones de inversion de spin, 
y acoplamiento spin-orbita [37,38]. Ademas estudiamos problemas que surgen al realizar 
la bosonizacion de teon'as de materia condensada en el marco de la integral funcional. El 
plan es el siguiente: en el capitulo |2l presentamos la bosonizacion en el marco operacional 
en forma detallada, y la aplicamos al estudio de modelos concretos de teorfas de muchos 
cuerpos en baja dimension. Para esto seguimos la bibliograffa estandar [17,39,40]. En el 
capitulo El analizamos las ambigiiedades que se presentan en la bosonizacion debido a los 
necesarios mecanismos de regularizacion que deben implementarse, particularmente cuando 
la teoria bajo estudio no posee la invarianza de Lorentz de las teorfas de campos usuales, y 
allanamos el camino para estudiar teorfas mas complejas de materia condensada mediante 
bosonizacion funcional [41]. En el capitulo IH atacamos uno de esos modelos: el modelo de 
Thirring no local con dos especies de fermiones, que es una version de teorfa de campos 
del modelo de Tomonaga-Luttinger con spin. Consideramos ademas el efecto de anadir al 
Lagrangiano terminos de inversion de spin. Obtenemos una accion bosonica efectiva que 
representa oscilaciones de densidad de carga y de spin de forma independiente; la funcion 
de particion resulta factorizada, dando lugar a la separacion spin-carga [42]. En el capftulo 
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El examinamos la aproximacion armonica autoconsistente y su utilizacion en el marco de la 
integral funcional y las teorias de materia condensada. Hallamos una formula para el gap 
de las excitaciones del sector de spin del modelo estudiado en el capitulo anterior, como 
funcion de potenciales arbitrarios de interacciones electron-electron de tipo dispersion hacia 
adelante [42]. En este capitulo proponemos ademas un nuevo metodo para determinar el 
parametro incognita asociado a la aproximacion y lo aplicamos al estudio del regimen de 
escala del modelo de Ising en 2D fuera del punto critico y en presencia de un campo magnetico 
h [43] . En el capitulo El calculamos funciones de correlacion en sistemas unidimensionales 
de electrones en interaccion en los que los grados de libertad de carga y spin se encuentran 
acoplados a traves de la interaccion spin-orbita. Este acoplamiento esta representado por 
una asimetria en el espectro libre de los electrones. Estudiamos fluctuaciones de tipo ondas 
de densidad de carga y spin, y de tipo superconductor singulete y triplete. Mostramos 
que la interaccion spin-orbita modifica los exponentes del decaimiento de las funciones de 
correlacion y el diagrama de fases del sistema. Ademas encontramos que susceptibilidades 
que eran finitas a bajas temperaturas, se vuelven divergentes cuando la interaccion spin- 
orbita es suficientemente intensa [44]. Para concluir, en el capitulo [71 reunimos los resultados 
mas destacados y las conclusiones. 

Finalmente queremos mencionar que las investigaciones realizadas en esta tesis se com- 
plementan con los estudios del modelo de Tomonaga Luttinger sin spin realizados en las Ref. 
45 y 46. En la primera se estudio el efecto de interacciones de tipo dispersion hacia adelante 
y umklapp no locales, y la generalizacion de la ecuacion del gap a estos casos, y en la segunda 
se considero la asimetria en el espectro libre, pero en el caso simplificado de electrones sin 
spin. 
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Capitulo 2 



Bosonizacion 



Presentaremos una derivacion exhaustiva de la tecnica de bosonizacion en el 
marco operacional. Seguiremos las referencias usuales para mostrar la equiva- 
lencia entre operadores fermionicos y operadores bosonicos. Este es el punto de 
partida para el estudio de teorias de materia condensada en una dimension. 

2.1 Equivalencia entre operadores fermionicos y boso- 
nicos 

2.1.1 Campos fermionicos 

Tomemos una teon'a que puede formularse en terminos de un conjunto de operadores de 
creacion y de aniquilacion fermionicos en una dimension espacial, que satisfacen relaciones 
canonicas de anticonmutacion 



Estos operadores estan etiquetados por los indices r, que distingue particulas que se mueven 
a la derecha (r = +1) o a la izquierda (r = —1), s que en general puede referirse a M 
especies de fermiones, por ejemplo en problemas de multiples cadenas, pero que usualmcnte 
se utilizara para indicar el spin electronico (s = +1 para spin para arriba y s = — 1 para 
spin para abajo), y un mdice discreto y no acotado k que denota el momento (o niimero de 
onda), de la forma 



Aqui L es la longitud asociada al tamaiio del sistema y 6b es un parametro que determina las 
condiciones de contorno del problema. k usualmente etiqueta las autoenergias del sistema 
libre (con eo correspondiente a la energia de Fermi ep). Que este mdice sea discreto y no aco- 
tado es un requisito indispensable para realizar una derivacion rigurosa de las identidades de 




(2.1) 





(2.2) 
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bosonizacion. Estas identidades son independientes de un problema especifico como puede 
serlo el modelo de Tomonaga-Luttinger, o el problema de Kondo; y de la relacion de disper- 
sion efc. Esto es posible porque diclias identidades son igualdades entre operadores, es decir, 
validas cuando actiian sobre cualquier estado del espacio de Fock. Son independientes en- 
tonces del Hamiltoniano, cuya forma detallada solo se vuelve relevante al calcular funciones 
de correlacion. Su aplicacion a modelos mas concretos sera analizada mas adelante. En esta 
seccion seguiremos en detalle la exposicion hecha en la Ref. 39 con dos diferencias mfnimas: 
un cambio en la normalizacion, para adecuarla a las aplicaciones de materia condensada mas 
usuales, y la inclusion explfcita del indice r, debido a que ciertos conmutadores dependen 
explicitamente de el (en la mencionada referenda se lo incluyo junto con el mdice s en un 
unico indice rj = 1, M). 

Comenzando por el dado conjunto de operadores de destruccion c^rs con las propiedades 
fj2.1|) y ()2.2j) definimos un conjunto de campos fermionicos de la siguiente manera: 

^ OO _ CO 

A,{x) = E ^-(^) = 7l E ^"^'^'^4., (2.3) 

fc=— oo k=—oo 

donde x G [— oo, oo] es la variable espacial. Sus inversas son 



Ck 



rs 



Los operadores iprs satisfacen las condiciones de contorno 



dxe-^'-'^'^rsix), cl^ = — dxe'^-'^iPlix). (2.4) 







Asix + L) = e'^'^ij^x) y ijl{x + L) = e-'^'^ijUx), (2.5) 

periodicas para 6b = j antiperiodicas para 6b = 1. Las ecuaciones 1)2.11) y ()2.2|) . junto con 
la identidad 

E e^"^ = 27r E 5(?/ - 27rm) (2.6) 
implican de inmediato las relaciones de anticonmutacion 

{lljrs{x),^l,Ax')} =Srr'6ss' ' ^' ' (2-7) 

^rs{x),^r'Ax)] = Ul{x),^l,Ax)] = 0- (2.8) 



Para x, x' G [0, L] o L ^ oo, y condiciones de contorno periodicas se reducen a las relaciones 
usuales para campos fermionicos. 

El vacio fermionico |0)o (llamado a veces mar de Fermi) se define en la forma 
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Cfcrs|0)o = para k>0{n>0) (2.9) 

cL|0)o = para < (n < 0), (2.10) 

es decir que es un estado tal que para todos los valores de r y s los niveles llenos mas altos 
corresponden a = y los vacios ma bajos a = 1. Respecto a este vacio se define la 
operacion de orden normal del producto de operadores ABC . . . como 

: ABC • • • : = ABC ... - o{0\ABC . . . |0)o, (2.11) 
para A, B, C, . . . G {c^rs] cj.^^}. Esta definicion es equivalente a agrupar todos los operadores 

Ckrs con A; > y todos los c^^^ con /c < a la derecha de los demas. 
El operador numero de particulas de tipo rs se define como 



N 



J2 ■^s(^krs-= Yl [4rsCfcrs -o(0|4,,Cfe,,|0)o • (2.12) 
fe=— oo fe=— oo 

Mediante la aplicacion de operadores de creacion y aniquilacion sobre el vacio construimos 
estados con distintas configuraciones de particulas y agujeros, |N), autoestados de los oper- 
adores Nrs, con autovalores N^s 

Nrs\^) = Nrs\^). (2.13) 

Designamos con la letra N al conjunto de los autovalores N^s para los diferentes r,s, y 
por abuso de lenguaje diremos que un estado de N particulas es un estado en el que hay 
Nrs particulas de tipo rs. Notese que es posible aniquilar particulas con k < (ya que 
justamente el mar de Fermi esta lleno hasta el nivel k = 0). Alternativamente en estc caso 
decimos que creamos un agujero con impulso k. Esto disminuye el autovalor N^s, que puede 
tomar asi valores negativos. 

El conjunto de autoestados con un dado N conforman el espacio de Hilbert de N particulas 
Ti[^. El espacio dc Fock JF se define como suma directa de los espacios de Hilbert con numero 
fijo dc particulas JF = XleN '^n- 

Entre todos los estados con el mismo N hay uno que posee menor energia, es aquel que 
esta lleno hasta un determinado nivel, y vacio de alii en mas. Este es el estado fundamental 
de ?^N, |N)o. Podemos dar una definicion mas precisa de este estado: 

|N)o = n^-lO)o' (2.14) 



donde 



'^^Nrsrs^{Nrs-^)rs ■ ■ ■ '^Ir-s para A^^s > 0, 

C^;' = { 1 paraA^,, = 0, (2.15) 

,C{Nrs+l)rsC{Nrs+2)rs ■ ■ ■ COrs para A^s < 0. 
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2.1.2 Campos bosonicos 

A partir del estado |N)o pueden construirse el resto de las excitaciones de N particulas. 
Definimos los siguientes operadores de creacion y aniquilacion bosonicos que cumplen diclia 
tarea, 



bi 



qrs 



I 

fin 



E 



k+q rs^krsi 



-'qrs 



E 



''k—q rs 



C-kr 



(2.16) 



donde G Z 



es un entero positive, y q = 27mq/L > 0. Estos operadores, al actuar sobre 
cualquier estado |N) crean una combinacion de excitaciones de parti'cula-agujero sobre ese 
estado con q unidades de momento mas (o menos), pero sin salirse de TYn- En este sentido 
son operadores que aumentan y disminuyen el momento. Su normalizacion se eligio de modo 
que satisfagan relaciones de conmutacion bosonicas 



bqrsi bq'r's' 



b^ b"^ 

"qrsy "q'r's' 



0, 



Nr. 



b^ 

'bq'r's' 



0, 



(2.17) 



bqrs, ^\ir's' 



Srr'^ss' I y^k^q'-qrs^krs C^+g' r^Cfc+g r-s 

V ^9^9 fe=-oo 



brr'^ss'^qq' | [' '^krs'^^rs ■ ■ Cfc_|^^ j-^Cfc+g^s • 

|0)o 



(2.18) 



Las ecuaciones ()2.17j) se pueden verificar facilmente, pero la derivacion de ()2.18|1 requiere 
cierto cuidado, como notaron por primera vez Mattis y Lieb [16]: para q ^ q' los dos 
terminos en la primera Ifnea ya estan ordenados normalmente (esto es porque sus valores 
medios de vacio son nulos) y pueden restarse trivialmente mediante un cambio k —>■ k — q' en 
el segundo termino, dando cero como resultado. Sin embargo, para q = q' antes de hacer la 
sustraccion debemos construir expresiones ordenadas normalmente, de otro modo estan'amos 
restando expresiones infinitas de un modo no controlado. Los terminos en la segunda Ifnea 
se cancelan, reemplazando en el segundo termino k ^ k — q (esto ahora si se puede hacer 
porque estan ordenados normalmente). La definicion del vacio (Ecs. ()2.9p y (j2.10j) ) implica 
que la diferencia en los valores de expectacion de la tercera Imea arroja como resultado 



nA ^ ^ 



— nn = I. 



(2.19) 



n.fe=— oo nj;=— oo 
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Notese que la construccion de los operadores bgrs ()2.16p y la derivacion de los conmu- 
t adores ()2.18|) descansa fuertemente en el hecho de que el conjunto de ks es infinito y no 
acotado por debajo. 

Es facil verificar que dentro de TYn, |N)o actiia como estado fundamental para las excita- 
ciones bosonicas: 

^grs|N)o = 0, para todo g, r, s. (2.20) 

Intuitivamente esto es claro: si |N)o es el estado fundamental entre todos aquellos que 
contienen N particulas, entonces no se le pueden quitar unidades de momento sin quitar 
parti'culas, es decir, sin salir de TCn- 

Es obvio que los estados excitados |N) que conforman el espacio de Hilbert de N particulas 
se pueden obtener actuando sobre |N)o con alguna funcion de los operadores fermionicos: 
|N) = /(c|,^5, CfcVs)|N)o. Haldane [18] mostro que tambien existe una representacion en 
terminos de los 6^^^. Mas especificamente, mostro lo siguiente: 

Teorema 1. Para cualquier estado |N), existe una funcion f{h'^) tal que 

|N) = /(6t)|N)o. (2.21) 

Esta es una afirmacion para nada trivial ya que los operadores 6^ crean complejas combina- 
ciones de excitaciones particula-agujero; y constituye el corazon de la bosonizacion debido a 
que implica una igualdad entre espacios de Fock bosonicos y fermionicos. Omitiremos aqui 
la demostracion, y remitiremos al lector a la mencionada referenda, y tambien a la Ref. 39. 

El estado fundamental |N)o sirve para definir una operacion de orden normal bosonica 
de un producto de operadores de tipo hq^s J ^Jrs manera analoga al orden normal fer- 
mionico ()2.11|1 . Mas aun, ambos son equivalentes, es decir que si un producto de operadores 
bosonicos esta ordenado normalmente de acuerdo al orden bosonico, entonces tambien lo esta 
de acuerdo al orden fermionico, y viceversa. Por este motivo se utiliza la misma notacion 
para ambos. 

Con los operadores bosonicos definidos en la Ec ()2.16|) podemos definir campos bosonicos: 



- ^e^^'^Vse"""/', dsi^) = -J2 -^e-'^'^rse'^'^^ (2.22) 

q>0 V ' 

y su combinacion liermftica 



(?>0 V 9 g>0 V "3 



{X) = ^rs{x) ^ = - 

— ^/^a 
g>0 V 1 



J2 -1^ [e'^'i-hqrs + e-'^'^^fe;,,) e-^'i'\ (2.23) 



Aqui a > es un parametro infinitesimal que regulariza divergencias ultravioletas que ocur- 
ren en ciertas expresiones y conmutadores no ordenados normalmente. Usualmente se toma 
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del orden del espaciado de red a ~ 1/k-p. Los campos asi definidos satisfacen las relaciones 
de conmutacion 



q[ir{x—x')~a] 



q>0 

= -5,,/5,,,,ln |l -e^t*"(^-^')-'^: 

L— >oo 



—Srr'^ss' In 



— [a — ir{x — x')] 



(2.24) 
(2.25) 

(2.26) 

(2.27) 



La Ec. ()2.26j) se obtuvo utilizando la expansion en serie de log(l — y). Aqui se ve claramente 
que a actiia como cut-off de la divergencia ultravioleta para x = x'. Estos conmutadores 
son utiles en la evaluacion del producto de operadores de vertice (exponenciales de campos 
bosonicos). Utilizando la identidad 



(2.28) 



para operadores A j B que conmutan con [A, B], obtenemos 



T \ 1/2 

iiptsix) iiprsix) _ i(ipts+<fr3)ix) [iiplsix),iiprsix)]/2 _ I \ J^rs(x) ('2 29) 

V27ray ^ ■ ' 

_ g-i(¥'rs+</'L){xOg[~i</'rs(a.'). -«'pL(2.-)]/2 _ i ^-~i4>rs(x) (2 30) 

Notese que estas formulas son validas para cualquier valor de L siempre que a sea suficien- 
temente chico (esto es asi porque para x = x' el limite L ^ cxd en fl2.26|) es equivalente a 
a — s> 0). Resulta interesante tambien la evaluacion del conmutador del campo 4>rs{x) con su 
derivada: 



A partir de aqui podemos obtener dos expresiones diferentes de acuerdo a como se tomen 
los limites para L ^ oo y a infinitesimal. Si queremos una expresion no periodica, para L 
grande, es conveniente hacer la suma geometrica, y posteriormente tomar los Ifmites dejando 
el limite a — >■ para el final: 



■iir- 



2tx 

T 



E 



2tt 

e L 



[ir{x—x')—a]nq 



+ 



2tt r 

e ^ ' 



-ir{x—x')—a\nq 



(2.31) 
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^°^> —Srv'Sss'^^Tcir- > —5rr'Sss''^'n'ir5(x — x'). (2.32) 

[x — x'Y + 

Notese que para tomar correctamente el limite L ^ oo en la primera Imea de las expresiones 
precedentes, se deben desarroUar los exponenciales hasta orden cuadratico en 1/L. Para L 
finite, en cambio, tomamos primero el li'mite a — en 1)2.311) . y utilizamos la identidad ()2.(i)) : 

(t)rs{.x),d^i(l)r's'{,x') 

donde el termino 1/L en esta ultima ecuacion aparece debido a la ausencia del termino 
= (g = 0) en la Ec. ()2.3H) . Finalmente podemos calcular el conmutador del campo (prs 
con si mismo, obteniendo 



-5rr'5ss'27rzr 



5{x — x' — nL) 



1 
L 



(2.33) 



(prs{x), (pr's'ix') °> 27iir5rr'Sss'G{x — x') doude e(a;) = < ' ^ ' (2.34) 

-I 10 si X = 0. 

2.1.3 Factores de Klein 

Los operadores h y Ij^ crean excitaciones dentro del espacio de Hilbert de N partfculas. 
Debemos definir entonces operadores que conecten espacios de Hilbert con diferente numero 
de partfculas, es decir, operadores escalera que aumenten o disminuyan el numero fermionico 
total, cosa que no pueden hacer los operadores bosonicos. 

Definimos los factores de Klein F^s y F^.^ como operadores con las siguientes propiedades: 
i) conmutan con todos los operadores bosonicos: 



h F^ 



^qrs 7 -^r's' 



F^ 

qrsi r's' 



para todo g, r, r', s, s', (2.35) 



y ii) su accion sobre un estado |N)o, es la de agregar una particula en el nivel mas bajo 
posible, y la de quitar una en el mas alto respectivamente: 



(Nrs+l)rs 



N)o, 



i^rs|N)o = CAr^,rs|N)o. 



(2.36) 
(2.37) 
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de este modo queda definida su accion sobre cualquier estado |N). En efecto, el estado |N) 
puede descomponerse de acuerdo a la Ec. ()2.21|) . y por lo tanto, 



i^,,|N) = /(6t)c^^^,,|N)o. 



(2.38) 
(2.39) 



Es decir, que el estado -Fj^g|N) (o Frs|N)) posee el mismo conjunto de excitaciones bosonicas 
que el estado |N), pero creadas sobre un estado con una particula mas (o menos). Asi 
definidos, los factores de Klein poseen las siguientes propiedades: 



FrsFlg = Fj^Frs = 1 (unitariedad); 

Fr's' I = 26rr'Sss' paia todo r, r', s, s'; 
I F^ , F^^, } = { Fr's' } = para r ^ r', s ^ s'; 



rs) ^ r's' 



5rr'Sss'Fjg, 



Srr' ^As' F'fa . 



(2.40) 
(2.41) 

(2.42) 

(2.43) 



Para probar la unitariedad es fundamental que el espectro del operador Nrs sea no acotado. 



2.1.4 Identidades de bosonizacion 

Con todas las definiciones y propiedades estudiadas estamos en condiciones de establecer 
igualdades entre operadores de campos bosonicos y fermionicos. La primera de ellas, la mas 
simple de derivar, establece una igualdad entre la densidad electronica ordenada normal- 
mente, y la derivada del campo bosonico dx(prs(yX): 



Prs{x) ^ : ^Ux)iJrs{x) := - J] e-'^- ■ ^UrsCkrs : (2.44) 

q k 

= ^ E (^^^'^ V. - ^~''''%rs) + \Y.-- ^Lc... : (2.45) 

q>Q k 

= - r-^dx(f)rs{x) + ^Nrs (para a 0). (2.46) 

ZTC L 

La segunda, relaciona el campo fermionico con el operador de vertice bosonico. Para 
derivarla debemos mostrar previamente la siguiente propiedad: 

Prop. 1. il)rs{x)\^)o es un estado coherente bosonieo 



18 



Mostraremos que dicho estado es un autoestado de bg^s Y por lo tanto posee una repre- 
sentacion como estado coherente. Para ello basta con calcular los conmutadores de b j 
con ip: 



bqr's'-,i>rs{x) = 6rr'5ss'0:qrix)'lprsix) 



(2.47) 
(2.48) 



donde aqr{x) = Estos conmutadores y la ecuacion ()2.2()j) implican inmediatamente 

que 

bgr's'1prs{x)\^)o = 6rr'Sss'aqr{x)lprs{x)\N)o. (2.49) 

Y por lo tanto, este estado posee una representacion como estado coherente bosonico [47]: 



iprs{x)\N)o = exp 



E 

.q>o 



0^qrix)b^qj.g 



ErsXrs{xm)o = C"^^- F,, A,, (x) | N)o 



(2.50) 



Aqui utilizamos la definicion del campo ()2.22|1 en la segunda igualdad. Hemos agregado 
el operador de fase A^^ que derivaremos en lo sucesivo; y el factor de Klein, que es necesario 
porque iprs remueve una particula del estado |N)o, cosa que los campos bosonicos b'^ no 
pueden hacer. Para obtener el operador A calculamos el siguiente valor medio de dos formas 
diferentes: por un lado. 



o{N\Fj^^rs{x)\M)o = o(N|i,(x)|N)o = Xrs{x) 



(2.51) 



donde hemos pasado adelante el factor de Klein F^s en ()2.50|) . ya que segiin su definicion 
()2.35|) conmuta con todos los 6^; utilizamos la unitariedad de los F's, y expandimos en serie 
el exponencial, quedandonos con el termino de orden 0, ya que o(N|6j^g = 0. 

Por otro lado, insertamos la descomposicion de Fourier ()2.3j) para ip{x) y la definicion 
del factor de Klein ()2.Hfi|) . y nos quedamos solo con el termino Uk = Nrs (o bien k = 



(N|F;>,,(x)|N)c 



Concluimos entonces que el operador Xrs{x) esta dado por 



\rs{x) 



(2.52) 



(2.53) 
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Para derivar las identidades de bosonizacion debemos estudiar la accion del campo iprs{x) 
sobre un estado |N) arbitrario (que segiin ()2.2H) puede escribirse como |N) = /({&gr's'})l'^)o)- 
Para ello utilizaremos las siguientes identidades [39], 

= fiibl's' - Srr'6ss'a;,ix)})^Prsix), (2.54) 

fiibU' - Srr'Sss'al^ix)}) = e-^"(^V({&;.J)e^^-(^\ (2.55) 
que se pueden mostrar facilmente a partir de la formula de Baker-Hausdorff 

e-^Ae^ = A+[A,B] + ^[[A, B]B] + . . . , (2.56) 

expandiendo en serie de Taylor la funcion / y empleando los conmutadores ()2.47|) y ()2.48|) . 
Podemos evaluar entonces iprs{x)\N) conmutando iprs{x) con /({fe^j.^}), insertando la repre- 
sentacion ()2.5()|1 y reordenando los factores: 

ijrs{xm)=^rs{x)f{{bUm)0 





- 6rr'Sss'a*qri^)})As{x)\N)o 


[por la Ec. 




= fiibl's' 


-6rr'6ss'al,{x)})e-'^'^^^^^Frs\rs{x)\N)o 


[por la Ec. 




-Frs^rs^X 


)e-^^'-(^^f{{bl,, - 5rr'5ss'al.{x)}mo 


[por la Ec. 




Frs^rs^X 


)^-ids{x) [e-'^-("V({&lrs})e'^^^^"T |N)o 


[por la Ec. 




-Frs^rsix 


)g-ic,L(xOg-i^.4-)/({5t^j)|N)o 


[por la Ec. 




Frs^rsi^X 




[por la Ec. 





Dado que | N) es arbitrario, y que todo estado del espacio de Fock es de esta forma, concluimos 
que las siguientes formulas de bosonizacion valen como identidades entre operadores en el 
espacio de Fock, y para todo L: 

ijrsix) = F,,A,3(a;)e-'^^=(")e-'^'-»(") (2.58) 
= i7;^_Le^'-¥(^--|^6)^e-^^^»(^)e-^'^-(^) [por la Ec. 1233] (2.59) 

= ^^i7;^e*"^(^--^^'')^e-^'^-(^'). [por la Ec. EM (2-60) 

\'2na 
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Por ultimo estudiaremos como representar un Hamiltoniano fermionico libre con una 
relacion de dispersion lineal. Mas especfficamente tomemos un Hamiltoniano de la forma 
(con h= 1) 

Ho = J2Hors (2.61) 

r,s 

con 

Hors = -irvrs / dx : V'j;^(a;)9^Vrs(a;) := ^ : cl^s'^krs ■ (2.62) 

La segunda forma se obtiene de la primera insertando el desarrollo de Fourier del campo 

La ecuacion ()2.21|1 implica que los b'^'s actuando sobre |N)o generan todo el espacio de 
Hilbert de N particulas. Esto significa entonces, que Hors debe tener una represent acion en 
terminos solamente de variables bosonicas. Para hallar esa represent acion, estudiemos el 
conmutador de fej^^ con Hq^s- 

= qblrs^rr'Sss'- (2.63) 

Ademas, dado que [-ffors, Nr's'] = para todo r, r', s, s', todo autoestado de N^s lo es tambien 
de Hors, en particular el estado fundamental de N particulas, |N)o. Su autovalor es 

i^ots = o(N|//ors|N)o = ^ (^^^ Nrs (iV,, + 1 - 5,) . (2.64) 

Se comprueba que la linica forma bosonica para ifors que reproduce las Ecs. ()2.63|) y 
es: 



Or 



Hors = J2 "l^rsKs + ^ ( Y ) + 1 - (2.65) 

<?>0 ^ ^ 

= -y[^'- (9^^rsix)f : +^ Mrs [Mrs + 1-6,). (2.66) 

Con esta ultima ecuacion completamos la derivacion de las identidades de bosonizacion, 
que valen para L finito. Para obtener expresiones con L — > oo basta con despreciar los 
terminos ~ 1/L. En este capitulo seguimos un enfoque constructivo, de modo que no es 
necesario verificar los conmutadores de los campos fermionicos o igualdad entre funciones 
de Green. A continuacion veremos las aplicaciones del proceso de bosonizacion, y como 
se vuelve extremadamente iltil para el estudio de complicadas teorias fermionicas en una 
dimension espacial. 
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ek 




2.2 Aplicaciones de la bosonizacion: el modelo g-ology 

Veremos a continuacion como aplicar la tecnica a una teoria de electrones interactuantes 
en una dimension. El punto de comienzo de las tcorfas de muchos cuerpos en materia 
condensada es un Hamiltoniano de la forma H — Hq-\- i^int, con 

^o=5]J]ek4A., (2.67) 

k s 

Hint /?''^'4+q.4-c.'4'.'4., (2.68) 

qkk' ss' 

donde V es el voliimen del sistema, y c|^^ y Cks son los operadores de creacion y aniquilacion 
de electrones, que satisfacen las relaciones de anticonmutacion canonicas 

|cks,Ckv} = ^kk'f^ss'- (2.69) 

Las cantidades f^^^'^' son los Uamados pardmetros de Landau, que describen la dispersion de 
dos particulas desdc un cstado inicial con numeros cuanticos (k, s) y (k', s) a un estado final 
con numeros cuanticos (k + q, s) y (k — q, s). En general k es el momento o pseudomomento 
de la particula, s su spin y q es el momento transferido. ek es la energia cinetica de los 
electrones medida desde el nivel de Fermi. 
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ek 



-n/a 




n/a 



k 



Figura 2.2: Relacion de dispersion linealizada alrededor de los puntos de Fermi i/cp- 



Como ya hemos comentado, en esta tesis estamos interesados principalmente en prob- 
lemas de electrones en una dimension espacial, y de aqui en mas nos restrigiremos a este 
caso. En la aproximacion de electrones casi libres, o de ligadura fuerte, la energi'a cinetica o 
relacion de dispersion se ilustra esquematicamente en la Fig. 12.11 En particular, en el mod- 
elo de Hubbard [48], el modelo realista mas simple que puede plantearse en una dimension, 
tenemos una expresion explicita para e^: 



donde t es la constante de intercambio (acoplamiento entre sitios vecinos) y a es el espaciado 
entre los atomos de la red. 

El espectro del modelo de Hubbard puede hallarse mediante la tecnica del ansatz de 
Bethe [49] (lo que a menudo en la literatura se denomina "resolver"el modelo). Sin embargo, 
esta tecnica no brinda resultados para las funciones de Green. Por este motivo se vuelve 
necesaria la aplicacion de otros metodos que permitan obtener una descripcion de los tipos 
de fluctuaciones que tienen lugar. Una de las tecnicas mas utilizadas para el calculo de 
funciones de correlacion es la bosonizacion. Como se ha mostrado en las secciones previas, 
en su formulacion mas usual (operacional) este procedimiento se basa en identidades entre los 
operadores fermionicos originales de la teon'a y ciertos operadores bosonicos que se pueden 
construir a partir de los primeros. En particular, el Hamiltoniano de los diferentes modelos 
resulta escrito completamente en terminos de fluctuaciones bosonicas. Sin embargo debemos 
hacer una aclaracion: el modelo de Hubbard originalmente se formula en la red, mientras 
que la bosonizacion se aplica a modelos continuos. Por ello, al estudiar este modelo, en 



efc = —2t cos ka 



(2.70) 
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algun punto debe tomarse dicho limite. Al hacerlo las propiedades que se derivan tienen 
validez para grandes distancias (ej. los decaimientos de las funciones de Green), o lo que es 
equivalente, para las excitaciones de baja energi'a. 

La superficie de Fermi en un metal estrictamente unidimensional consiste en dos puntos, 
+/cf y —kp] en su vecindad podemos linealizar la relacion de dispersion ()2.7()|1 [50]: 

ek = VFi\k\-h) (2.71) 

Esta aproximacion, en principio, es razonable en un rango finite alrededor de los puntos de 
Fermi. Sin embargo, los tratamientos matematicos se simplifican enormemente si tomamos 
esta version linealizada para todos los valores de k entre — oo y +00, es decir, si reemplazamos 
el espectro libre de la Fig. l2.1l por el de la Fig. 12.21 For otro lado, de acuerdo a lo dicho mas 
arriba, solo estamos interesados en excitaciones de baja energi'a, a las que contribuyen estados 
proximos a la superficie de Fermi; de modo que la inclusion de los estados adicionales por 
efecto de la linealizacion es despreciable en este regimen. Esto ultimo se verifica a posteriori 
al estudiar los efectos producidos por la curvatura de banda, es decir incluyendo terminos 
cuadraticos y cubicos en la relacion de dispesion. Se puede mostrar que las contribuciones 
de estos terminos son irrelevantes frente a las del termino lineal. 

La linealizacion genera dos ramas bien definidas en la relacion de dispersion. Los elec- 
trones que pertenecen a la rama que contiene al punto +kp (— A^f) se mueven hacia la 
derecha (izquierda), a los operadores que los representan los denotaremos c|^_^,p)^5 y C(^k~kp)Rs 

{c'(^-k-kp)Ls y C(-k-kp)Ls)- Este conjunto de operadores asi definidos satisface los requisitos de 
ser un conjunto infinito y no acotado, y los identificamos inmediatamente con los descriptos 
en las Ecs. (12.11) y ()2.2|) . En termino de estos operadores, el Hamiltoniano libre se escribe 



Ho = Yl vpk clj^s^kRs + c\LsCkLs : j = -IVY I ■■ y^Rsdx^Rs - ipis^xipLsj ■ ■ (2.72) 

ks s "'O 

Para obtener la primera igualdad cambiamos A; — en el segundo termino, y posterior- 
mente efectuamos la traslacion k —>■ k + kp. Para obtener la segunda, utilizamos la definicion 
de los operadores de campo fermionicos ()2.3|) . 

Los terminos de interaccion se pueden clasificar en cuatro tipos diferentes mostrados en 
la Fig. 12.31 Los electrones pertenecientes a ambas ramas se distinguen mediante Ifneas 
punteadas y solidas. El proceso con constante de acoplamiento gi corresponde a dispersion 
hacia atras, y posee una transferencia de momento de 2kp. Los procesos con constantes g2 
y g4 son de dispersion hacia adelante, su transferencia de momento es nula. Por ultimo el 
proceso con constante 5^3 es de tipo umklapp, y su transferencia de momento es de Akp. Este 
ultimo proceso solo es importante cuando estamos con un Uenado medio, es decir, pensando 
en una situacion del tipo del modelo de Hubbard, cuando tenemos un electron por sitio. En 
ese caso, ikp es igual al vector de red reciproco. 
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Figura 2.3: Procesos que intervienen en la dispesion del gas de electrones unidimensional: 
gi dispersion hacia atras; g2 y Qa dispersion hacia adelante; umklapp. 



Al introducir el spin electronico, para cada uno de estos procesos aparecen dos variantes, 
de acuerdo a la oricntacion relativa de los electrones incidentes. Si ambos poseen spines 
alineados, le agregamos un submdice || a la constante de acoplamiento; en cambio, si los 
spines son antiparalelos, agregamos el submdice _L. Notese que estos procesos no invierten 
el spin, es decir, el spin de las particular finales es identico al de las parti'culas iniciales. 
La generalizacon de este modelo al caso en que las interacciones pueden cambiar el estado 
de spin es uno de los objetivos de esta tesis (ver Capitulo 4). Volviendo al caso presente, 
la expresion para el Hamiltoniano de interaccion que describe estos procesos, en espacio de 
coordenadas es 
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ss' 

+ {92\\Sss' + 92±Ss-s') pLspRs' 

ss' 

^ rL 

i^^W^"'' ^ 33±Ss-s') i^-rsi^rsi^lrs'i^rs' 

rss' 

1 

+ {diW^ss' + g4±Ss-s') PrsPrs'- (2.73) 

rss' ^ 

Prs esta definido en la Ec. ()2.44p . La teoria definida por las dos ecuaciones anteriores, version 
continua del modelo de bajas energias de un ensemble de electrones, se conoce popularmente 
como modelo de "geologi'a" , o "g-ology" [50] . La bosonizacion del Hamiltoniano es inmediata 
aplicando las identidades ya mostradas en la seccion previa. Por simplicidad tomaremos 
L ^ oo, y 5(311 = 5(111 = 0, g3± = 53, gi± = gi. Para Hq usamos las ecs. (I2.(il|). (12.621) y (12.661): 

Ho = 'fj:l^(9.^rs)\ (2.74) 

rs 

mientras que para ifjnt utilizamos p.44j) - p.46j) para los terminos de dispersion hacia ade- 
lante, y (j2.6()j) para los de dispersion hacia atras y umklapp. El resultado es 



-«<PL,- 



s 

+ ^ I dx {g2\\Sss' + g2±Ss-s') -7^^d^(j)Lsdx4>Rs' 
ss' ^ 

+ I dxg3e'^--'e-''t''-'e''''-^'-'e-"^^'-' 

rs 

1 V f 1 

+7:^ dx {g^\Sss' + g4±Ss-s') j—r^dx4>rsdx4>rs'- (2.75) 

^ rss' J ^^'^> 

En los terminos de dispersion hacia atras y umklapp no tuvimos en cuenta los factores de 
Klein porque sus contribuciones pueden despreciarse en el limite L — »• 00 [17]. Introducimos 
a continuacion los campos Oy y (py con u = p,a. 

^ [9, -r(Pp + s {6, - r0,)] . (2.76) 
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Para entender el significado fisico de estos campos, estudiamos los operadores de densidad 
de carga y de spin, definidos del siguiente modo: 



P 



E 



Pr 



a 



E 



SPr 



(2.77) 



Utilizando la equivalencia de bosonizacion ()2.46p . podemos encontrar expresiones para estos 
operadores en terminos de los campos recientemente definidos: 



P 



'^d^(j)p; a = y ^<9^0<x- (2.78) 

Entonces los campos (/)p y estan relacionados con la densidad de carga y de spin respec- 
tivamente. Mas aun, utilizando los conmutadores del campo (prs con si mismo ()2.34|) y con 
su derivada ()2.3H) es posible mostrar que 0p y 9p, y sus derivadas conmutan con (pa Y Y 
sus derivadas. Ademas si definimos 11^ = 0^9^, liallamos que 



i6f,uS{x - x') 



(2.79) 



que constituyen las relaciones de conmutacion canonicas para campos bosonicos. Introduci- 
mos ademas las constantes 



92 
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las constantes de dureza 



(fi'2|| ± 92a) , 
(fl'4|| ± 5'4±) , 

^nvY + 91-92 
t^vy + 91 + 92 



Y las velocidades renormalizadas 



v/(vri;F + 9l - 92) {t^^f + 9l + 92. 



(2.80) 
(2.81) 

(2.82) 
(2.83) 



Con estas definiciones el Hamiltoniano total se escribe como H = Hp + H^, donde 



dx— 
2 



dx— 
2 



Kp{djpf + ^{d^ct>pf 



2^1 

, 2^3 
+ 77^ — ^cos 



(2.84) 



(2.85) 



De este modo el Hamiltoniano fermionico original queda escrito completamente en termino 
de dos campos bosonicos (pp y (pa que representan oscilaciones independientes de densidad 
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de carga y spin respectivamente, y que se propagan con velocidades diferentes. Ademas los 
campos de spin conmutan con los campos de carga, y lo mismo ocurre con los Hamiltoni- 
anos de cada sector. Estas caracten'sticas dan lugar a una propiedad fundamental de este 
tipo de modelos en una dimension: la separacion spin-carga. Por ultimo mcncionemos que 
las funciones de Green del sistema admiten una factorizacion semejante, en tcrminos de un 
factor de carga y un factor de spin. Uno de los objetivos de esta tesis es analizar de que 
manera este fenomeno emerge en un estudio de la teoria en el marco de la integral funcional, 
y como se manifiesta al nivel de las funciones de Green y del diagrama de fases la ruptura 
de la separacion spin carga mediante la presencia de interacciones spin-orbita (ver Gapitulo 
6). 
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Capitulo 3 



El determinante fermionico no 
covariante y su relacion con los 
liquidos de Luttinger 

En este capitulo consideramos el procedimiento de bosonizacion cn el marco de 
la integral funcional. Esto nos conduce naturalmente al estudio del determinante 
fermionico asociado a una Teoria Cuantica de Campos no covariante, utilizada 
para describir un sistema no relativista en (1 + 1) dimensiones. Explotando 
la libertad que brinda no mantener la invarianza de Lorentz, determinamos el 
operador regulador correspondiente al metodo del nucleo del calor (heat-kernel) 
que permite reproducir la relacion de dispersion de las excitaciones bosonicas y 
los exponentes criticos correctos del modelo de Tomonaga-Luttinger. Ademas 
derivamos el Hamiltoniano del modelo bosonizado funcionalmente y las corre- 
spondientes corrientes. De este modo establecemos la regularizacion precisa 
mediante heat-kernel, que conduce a un complcto acuerdo cntre el abordaje op- 
eracional a la bosonizacion de modelos de materia condensada, y su alternativa 
mediante integrales funcionales. Estos resultados son parte de las contribuciones 
originales a esta tesis [41]. 

3.1 Introduccion 

Los determinantes fermionicos juegan un rol central en las formulaciones modernas de las 
Teorias Cuanticas de Campos (QFT's). Como es bien sabido, emergen naturalmente al con- 
siderar funcionales generatrices asociadas a campos fermionicos en el marco de la integral 
funcional [51]. En los ultimos veinte anos ha sido especialmente fructifero el estudio de de- 
terminantes fermionicos en (1 + 1) dimensiones. La observacion de Fujikawa concerniente 
a la no trivialidad del jacobiano asociado a las transformaciones quirales en las variables 
fermionicas [30,52,53] fue aplicada al caso en (1 + 1) dimensiones, y condujo a avances 
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significativos en nuestra comprension de los Uamados modelos de juguete paradigmaticos, 
tales como la electrodinamica cuantica en dos dimensiones (QED2), el modelo de Thirring, 
y sus versiones no abelianas [31,54,55]. De hecho, basado en un tratamiento adecuado 
del determinante fermionico, se desarroUo una tecnica de bosonizacion mediante integrales 
funcionales [56]. El punto crucial se encuentra en el calculo del mencionado jacobiano. 
Un calculo ingenuo arroja un resultado mal definido y se vuelve necesario implementar un 
procedimiento de regularizacion. En teorias de gauge con fermiones de Dirac, es natural 
considerar un esquema de regularizacion que preserve la invarianza de gauge. Por otro lado, 
cuando los campos vectoriales que se encuentran presentes en la teoria son solo campos aux- 
iliares (usualmente introducidos a traves de una transformacion de Hubbard-Stratonovich), 
se puede elegir un regulador mas general [57,58]. El modelo de Thirring [27] y el mod- 
elo de Schwinger quiral [59, 60] constituyen ejemplos en los que tienen lugar este tipo de 
ambigiiedades en la regularizacion. 

El tema de la regularizacion del jacobiano de Fujikawa, su relacion con contraterminos 
locales, y su rol en el analisis de anomalias cuanticas ha sido extensivamente examinado 
en la literatura [61]. En todos los casos, los modelos bajo estudio son QFT's relativistas, 
es decir, teorias covariantes de Lorentz. Sin embargo, en ciertas situaciones relevantes, el 
interes recae en teorias de campos no covariantes. Es el caso del analisis de sistemas de 
electrones unidimensionales que pueden estudiarse mediante el modelo g-ology. En este 
contexto, la bosonizacion funcional, alternativa al enfoque operacional usual, fue sugerida 
por primera vez por Fogebdy [33] y sucesivamente elaborada por Lee y Chen [34]. La 
conexion explfcita entre la bosonizacion funcional que condujo a una accion efectiva para la 
dinamica de las excitaciones bosonicas colectivas y el jacobiano de Fujikawa fue establecida 
por primera vez en la Ref. 35. Pero aun en este caso se empleo una regularizacion covariante, 
tomada de la teori'a de campos relativista. Como resultado las expresiones generales para 
las relaciones de dispersion de los modos bosonicos y los exponentes de decaimiento de las 
funciones de correlacion en terminos de las constantes de acoplamiento iniciales del modelo 
fermionico no acordaron con las obtenidas mediante bosonizacion operacional usual. En este 
capftulo mostramos que el origen de este desacuerdo se encuentra en el tipo de regularizacion 
escogido para calcular el jacobiano de Fujikawa. Dado que se viola la invarianza de Lorentz 
subyacente, parecen posibles un niimero arbitrario de esquemas de regularizacion. Solo uno 
de ellos conduce al resultado usual para los modelos g-ology y Tomonaga-Luttinger. En este 
punto quisieramos enfatizar que no fuimos capaces de hallar un principio fisico que sirviese 
de guia para elegir a priori entre diferentes esquemas de regularizacion, al estilo del principio 
de preservacion a nivel cuantico de simetrias que se hallan en la teon'a bajo estudio a nivel 
clasico, como lo son las ya mencionadas simetrias de gauge o de Lorentz. Sin embargo, hasta 
donde sabemos, tal principio tampoco ha sido identificado en el marco de la bosonizacion 
operacional de teorias de materia condensada. Por supuesto, este es un aspecto importante 
que merece futuras investigaciones. 

El plan de este capitulo es el siguiente. En la seccion 13.21 presentamos el modelo y 
expresamos su funcional generatriz en terminos de un determinante fermionico. En la seccion 
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13.31 con el objeto de clarificar la discusion comenzamos con un trazado de los pasos principales 
del enfoque desacoplante a la bosonizacion, y los resultados que se obtienen al emplear 
una regularizacion estandar invariante de Lorentz. Aqui incluimos dos subsecciones donde 
presentamos dos tipos de regularizacion diferentes: el metodo point-splitting, y el metodo 
heat-kernel. En este ultimo caso determinamos la forma precisa del operador necesario para 
obtener la respuesta correcta para las relaciones de dispersion y los exponentes. En la seccion 
13.41 mostramos como derivar, en nuestro marco de bosonizacion funcional, el Hamiltoniano 
bosonico y las correspondientes corrientes bosonizadas. Finalmente discutimos brevemente 
la conservacion de la corriente. En la seccion ITHl reunimos los resultados y las conclusiones. 

3.2 El modelo y el determinante fermionico 

Consideraremos una version no covariante del modelo de Thirring definido por el Lagrangiano 
euclideo 

L = ^i^^ - — (3.1) 

donde Vq and Vi son las constantes de acoplamiento y las derivadas estan redefinidas para 
incluir a la velocidad de Fermi: 

d 

di =t'F7^- (3.3) 

OXi 

Notese que vp juega el rol de la velocidad de la luz en QFT que usualmente se toma unitaria. 
Para f f = 1 y Vq = Vi = 1 tenemos el modelo de Thirring usual (la constante se 
incluye para facilitar la comparacion con los resultados invariantes de Lorentz). La corriente 
fermionica se define como 

= ^7^-0, (3.4) 
la cual satisface la ley de conservacion clasica 

0^3, = 0. (3.5) 

La funcional generatriz es 



Z[S] = j VipViP exp - j d^x{L+j^S, 



(3.6) 
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Por medio de una transformacion de Hubbard-Stratonovich, puede ser escrita en la forma 



[S]^J^j VA,detp[A] exp d'x £yV^-^{x - y)(gA^- S,){x)(gA^ - S^)(y) 



donde 

y 



Ip [A] =i^ + g4, 



det p[A]= I VipVi! exp 



- / d^x2pp[A]ij 



(3.7) 



(3.8) 



(3.9) 



3.3 Enfoque desacoplante de la bosonizacion 

Habiendo expresado la funcional generatriz en terminos de un determinante fermionico, 
trazaremos ahora un esquema del metodo desacoplante, que se encuentra en la base del 
enfoque funcional de la bosonizacion [31,54-56]. En (1 + 1) dimcnsiones espacio-tcmporalcs, 
el campo vectorial A^^ puede descomponerse en sus partes transversal y longitudinal del 
siguiente modo: 

A^^-(l/g){e^,d,(l>-d^ri), (3.10) 

donde r] (0) es un campo escalar (pseudoescalar) . Notemos que si realizamos la siguiente 
transformacion en los campos fermionicos 



-Qth5<l>-iri]-^ 



(3.11) 
(3.12) 



con t un parametro real, entonces la densidad lagrangiana fermionica cambia como 

iPp[A]ij = xPt[A]x (3.13) 



donde 



P[A]=p[{l-t)A]. 



(3.14) 



Como fue observado por primera vez por Fujikawa [30], el jacobiano asociado al men- 
cionado cambio en las variables fermionicas es no trivial, y depende de los campos ((> y rj: 



det(i^ + g4) = J[(j), rj; t] det(i^ + g{l - t)4). 



(3.15) 
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Debe notarse que para t = 1 los grades de libertad bosonicos y fermionicos se desacoplan 
completamente. Puede mostrarse que 



J[0, T]] 1] = J = exp 



con 



uj{t) dt 



(3.16) 



u{t) = -tTpt[A]-'g4 



lim tr 



D 



(fxIpt[A] \x,y)g4{x), 



(3.17) 



donde tr^ indica la operacion de traza en el espacio de Dirac. Todas estas formulas guardan 
una gran analogia con las correspondientes a una QFT covariante. Mas aun, la linica difer- 
encia entre ellas es la presencia de vp en lugar de la velocidad de la luz, aunque esto tiene 
consecuencias no triviales. La ultima ecuacion debe ser regularizada, de otro modo aparecen 
divergencias, como es obvio al tomar el limite y ^ x. En QFT, cualquier regularizacion 
aceptable tiene que ser invariante de Lorentz. En el presente caso, no tenemos esa limitacion 
por dos razones: i) la invarianza de Lorentz esta rota desde el principio dado que estamos 
en una teoria no relativista; ii) Hay una covarianza remanente: la teoria sin interacciones es 
invariante con respecto al grupo de Lorentz, donde la velocidad de la luz se ha reemplazado 
por la de Fermi, pero esta es una simetria artificial, y no existe razon para respetarla. Mas 
aiin, al tomar Vq ^ V\ {g2 7^ Qa) queda explicitamente rota. Antes de utilizar la libertad 
que proviene de la ausencia de covarianza, seria instructivo rever los resultados obtenidos 
previamente eligiendo una regularizacion invariante de Lorentz [35] . Utilizando un regulador 
de la forma 



{m]Pt[A]^ + Pt[A]^mA])/2, (3.18) 

que fue propuesto por primera vez por Fujikawa en su analisis de las anomalias covariante y 
consistente [61], se obtiene 

Jco. = exp j d'x [(9i0)2 + {do^f] I , (3.19) 

donde a es un parametro vinculado a posibles ambigiiedades en la regularizacion. Para a = 1 
se obtiene una regularizacion invariante de gauge. Aunque el modelo de Thirring no posee 
invarianza de gauge local, en el presente contexto estamos interesados fundamentalmente 
en la invarianza de Lorentz y podemos fijar a = 1 sin perder generalidad. Insertando el 
jacobiano anterior en la funcional generatriz, absorbiendo el determinante fermionico libre, 
lo que resulta del procedimiento de desacople en un factor de normalizacion, y expresando 
en terminos de (p j 1] segiin la Ec. ()3.10p . se obtiene una accion bosonizada. En el contexto 
de la materia condensada estos grados de libertad bosonicos se interpretan como campos 
asociados a oscilaciones de densidad de carga. De esta accion bosonica derivada a traves 
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cov 



de una regularizacion que preserva la covarianza se puede facilmente calcular la relacion de 
dispersion correspondiente 

Po + ^covP? = (3.20) 

donde 

^vl) r- 3.21 

('*+^) 

Debemos enfatizar aqui, que solo para Vi = esta velocidad, Uamada velocidad renormaliza- 
da, acuerda con el valor obtenido usando bosonizacion operational en materia condensada, 
que es 

^ - [^F--^j[vF + -^j. (3.22) 

A continuacion describiremos dos metodos diferentes para regularizar el jacobiano, que 
no preservan la invarianza de Lorentz y que permiten obtener la accion bosonica efectiva que 
conduce a la respuesta correcta para las relaciones de dispersion. 

3.3.1 Metodo Point-splitting 

Como es bien sabido, el metodo de regularizacion point-splitting rompe la invarianza de 
Lorentz explicitamente. Consiste en una prescripcion para tomar el Hmite y — > x antes 
mencionado al definir 

\im pt[A]-\x,y) = Mm + \im)pt[A]-\xo,x,;xo,Xi + e), (3.23) 

es decir, tomando un Hmite simetrico en la variable espacial. Necesitamos entonces la funcion 
de Green del operador de Dirac, que satisface 

p,[Aim]-\x,y)^5'ix-y) (3.24) 
Como es usual, proponemos el ansatz 

pt[A]-\x, y) = e(^-*)[^5^(^)+^^(^^l ^0(2:, y) e(i-*)[^5'^(2')-'''(2')], (3.25) 
donde Go es la funcion de Green del operador de Dirac libre: 

i^,Go{x,y) = 5\x-y). (3.26) 
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Con esta receta, hallamos el siguiente resultado para la ecuacion ()3.23|) : 



\iTcypt[A]-\x,y) = --^(1 -t)7i9i [750(x) - 2r/(x)] (3.27) 
y entonces, el Jacobiano (Ecs. ()3.1(jj) y ()3.17j) ) esta dado por 

J = exp J d^x [id,<PY - id,T]Y - 29iMr/] | . (3.28) 

La funcional de vacio puede entonces escribirse como 

Z[S = 0]=Af j V(t>Vr]e-^''' (3.29) 

donde Af es un factor de normalizacion que incluye el determinante fermionico libre (inde- 
pendiente de las interacciones) . Tambien hemos definido Sei, que en el espacio de momentos 
toma la forma 

Sef = I 7^ [mM{-p) + V{P)BV{-P) + 20(p)Cr7(-p)] . (3.30) 



con 



A _ 2 2 ( ^ 1\ pI 

^=^^^°^K2?T^-2?T^-2;^) 



(3.31) 
(3.32) 
(3.33) 



El contenido fisico del modelo puede extraerse de Sef que describe la dinamica de los 
modos colectivos del sistema. Cuando el modelo fermionico original esta relacionado al 
modelo de Tomonaga-Luttinger utilizado para el estudio de sistemas electronicos en una 
dimension [17,40,62,63], estas excitaciones colectivas corresponden a oscilaciones de densidad 
de carga (plasmones). Su relacion de dispersion puede obtenerse de los ceros del determinante 
de la matriz 



A C 
C B 



(3.34) 



El resultado es 



Po + ^'p? = (3.35) 
donde v es la velocidad renormalizada de los modos de densidad de carga dada por la Ec. 
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3.3.2 Metodo Heat-kernel 



Otro modo popular de tratar la regularizacion de los determinantes fermionicos es el metodo 
heat-kernel [30,64]. En este esquema J se regula insertando un operador de la forma e~^/*^ , 
R es un operador definido positivo, y M es un parametro que juega el rol de una masa, y 
que se deja fijo en los calculos intermedios. El li'mite — > oo se toma al final. De nuevo, 
debemos mencionar que en contextos de QFT estandar el operador R puede elegirse entre 
aquellos compatibles con la invarianza de Lorentz (dejemos de lado, por el momento, otras 
posibles simetrias), por ejemplo R = pt[A]'^. Aqui no tenemos esa limitacion, y nuestro 
objetivo es hallar la forma precisa de R que conduce a una accion efectiva que contenga las 
relaciones de dispersion deseadas. 

Comenzamos por reescribir la ecuacion ()3.17p como 



La operacion de traza esta mal definida, y necesita ser regularizada. Definimos nuestra 
u! regularizada como 



u;{t)n = lim tr {MA]-' [M ' ^v) PM] + PM] M + ^v)] e"^/*'') . (3.37) 



La eleccion de R siempre se encuentra dictada por consideraciones fisicas, por ejemplo, si 
consideramos una teoria de gauge, debemos tener en cuenta regularizaciones que no destruyan 
la invarianza de gauge a nivel cuantico. Esto es usualmente realizado tomando R = pt[A]'^, 
donde es un campo de gauge. Aqui el modelo bajo estudio no es una teoria de gauge 
y por lo tanto tenemos mas libertad de elegir el regulador. Emplearemos un operador de 
la forma R = ptlB]"^, donde i?^ es un cierto campo vectorial a ser determinado. Podemos 
escribir uj{t)ji como uj{t)ji = uJo{t) + uj^c(t) donde 



Aqui el subindice indica el termino que habriamos obtenido si hubieramos empleado la 
propiedad ciclica de la traza en la ecuacion ()3.36j) . El subindice nc se refiere a un termino 
"no ciclico" (esta cuestion se discute en detalle en la Ref. 58). Las expresiones finales para 
estos dos terminos son 




(3.36) 




(3.38) 



(3.39) 




(3.40) 



(3.41) 
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En este punto, un calculo directo nos permite comprobar que tomando 



5o =Ao (3.42) 
Bi = - Au (3.43) 

se arriba al mismo resultado obtenido en la subseccion previa Ecs. fl3.28|) - ()3.34|) ). De este 
modo, hemos hallado una forma explicita para el operador que regula el jacobiano asociado 
a un determinante fermionico no covariante. A su vez, esta forma conduce a las relaciones de 
dispersion correctas para el modelo de Tomonaga-Luttinger. La derivacion de este jacobiano, 
y del operador regulador era una de las principales motivaciones de este capftulo, y constituye 
uno de los aportes originales de esta tesis. 



3.4 El Hamiltoniano bosonizado y las corrientes 

Hasta este punto hemos trabajado en la formulacion lagrangiana. Sin embargo en aplica- 
ciones de materia condensada, el marco Hamiltoniano es usualmente el preferido. Es deseable 
entonces mostrar la derivacion del Hamiltoniano usual para sistemas electronicos unidimen- 
sionales en el marco de la bosonizacion funcional discutido en este capftulo, es decir, la forma 
bosonica del modelo de Tomonaga-Luttinger [17,40,62,63]. El otro punto que tratamos en 
esta seccion es la forma bosonica de las corrientes fermionicas originales (densidad de carga 
y corriente electrica), y sus leyes de conservacion. 

Teniendo en cuenta la expresion para det p [A] calculada en las secciones precedentes, y 
la relacion entre los campos (p j rj j e\ campo A^ (Ec. (l3.10|) ). podemos expresar la funcional 
generatriz ()3.7j) en terminos del campo A^: 



Z[S]=Af j PA^exp (^-^ j (fx(fyA^{x)D^,{x-y)A,{y)^ x 

exp (^-^ j d'^^d'^yS^,{x)V^~^{x-y)S^{y)^exp (^-^ j d'^^ d^V A^{x)V^~){x - y)S^{y) 

(3.44) 

donde D^^, esta dado en el espacio de Fourier por 

DM= ''''2) + (j ?V (3.45) 

El campo A^ se desacopla de la fuente S'^ a traves del procedimiento usual de efectuar la 
traslacion 

A,^A, + (3.46) 
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obteniendo 



Z[S] J VA^e^p(^~J d^xd^yA^{x)D^,{x-y)A,{y)^ x 

exp ^ J d^xd'^y S^{x)A-^{x - y)S^{y) 
con y la constante de dureza K dados por (en el espacio de Fourier) 



K 



, (3.47) 



(3.48) 



(3.49) 



VF + g'^Vo/n' 

Esta constante gobierna los exponentcs del decaimiento de las funciones de correlacion. A 
continuacion podemos multiplicar y dividir por 

J VA^e^p J d^xd^yA^{x)A^,{x - y)A,{y)^ , (3.50) 
y efectuar una traslacion adicional en el campo ^4^ 



(3.51) 



para obtener 



Z[S]=^^ J P^^expj^-- j d^xd^yA^{x)A^,{x-y)A,{y)+ / d'^x S^{x)A^{x] 
Finalmente, definiendo los campos (fi y 9 del siguiente modo: 



(3.52) 



A^ = ^dJ, 



Uegamos a la funcional generatriz 



(3.53) 
(3.54) 



Z[S] =M / V^VOei^p ( / dxdr ^{d^^f + vK{djf + 2idjdr^ 



X 



exp dxdr [—So^a;(f/^/^^ + ^Sl^x9/^/^^\^. (3.55) 
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Naturalmente identificamos entonces el campo 93 con el modo de densidad de carga del 
sistema y U — d^O como su campo canonico conjugado. Mas aiin, los primeros dos terminos 
en la accion cuadratica de la expresion previa pueden ser identificados con el Hamiltoniano 
del sistema: 



^-fdx [^{d^vf + vK(ajf] , (3.56) 



« 2 

el cual coincide exactamente con el Hamiltoniano obtenido utilizando bosonizacion opera- 
cional estandar [17,40,62,63]. Ahora, por derivacion funcional obtenemos la forma bosonica 
de las corrientes 



jo = —F^ dxV (3.57) 



Ji = ^n, (3.58) 

que son identicas, por supuesto, a las halladas en el enfoque operacional. Es importante 
enfatizar que estas corrientes no obedecen la ecuacion de continuidad. Siguiendo la Ref. [65], 
se introduce una corriente electrica fisica j, que es en general diferente de ji. La densidad de 
carga es identificada con j'o (jo — p)- La corriente fisica se determina al imponer la ecuacion 
de continuidad: 

Obtenemos 

j^^vKIi. (3.60) 
Notese que solo para Vi — Q {g^ — en el lenguaje de Tomonaga-Luttinger) se Uega a 



3.5 Conclusiones 

En cstc capitulo consideramos un determinante fermionico asociado a teorias cuanticas de 
campos no covariantes. En particular estudiamos el determinante que aparece al implemen- 
tar la bosonizacion mediante integrales funcionales basada en el desacople del determinante 
fermionico a traves de cambios apropiados en las variables de integracion. El modelo anal- 
izado (una version no covariante del modelo de Thirring) se ha utilizado previamente para 
dcscribir sistemas de electrones altamente correlacionados en una dimension (liquidos de 
Luttinger). 
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En el contexto del metodo de regularizacion del niicleo del calor, explotando la libertad 
originada en la no covarianza, determinamos el operador que actiia como regulador y que da 
lugar a una accion bosonica, a una relacion de dispersion, y a un decaimiento de funciones 
de correlacion correctos en terminos dc las constantes de acoplamiento, es decir, en completo 
acuerdo con aquellas que se obtienen, y son bien conocidas en el marco operacional. Calculos 
previos mediante integrales funcionales habian hecho uso de reguladores covariantes tornados 
de teorias de campos relativistas, y daban un espectro y exponentes correctos solo para 
valores particulares de las constantes de acoplamiento. 

Finalmente mostramos como derivar el Hamiltoniano bosonizado y las corrientes , que 
coinciden con los obtcnidos mediante bosonizacion operacional estandar. De estc modo 
fuimos capaces de establecer el regulador heat-kernel que brinda completo acuerdo entre los 
enfoques operacionales y mediante integrales funcionales de la bosonizacion del modelo de 
Tomonaga-Luttinger. 
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Capitulo 4 



Interacciones de inversion de spin en 
el modelo de Thirring no local 

Extendemos una version no local y no covariante del modelo de Thirring con 
el objeto de describir sistemas de muchos cuerpos que poseen interacciones de 
inversion dc spin. Introduciendo un modelo con dos especies de fermioncs cvita- 
mos el uso de bosonizacion no abeliana, necesario en un enfoque previa. Obten- 
emos una expresion bosonizada para la funcion de particion, que describe la 
dinamica de los modos colectivos del sistema. Los resultados de este capitulo 
constituyen un aporte original de esta tesis [42]. 

4.1 Introduccion 

En el capitulo anterior estudiamos el modelo de Tomonaga-Luttinger mediante bosonizacion 
funcional. Como ya se ha mencionado este es un modelo en el que se desprecia el spin 
electronico y que solo contempla interacciones de dispersion hacia adelante, que poseen 
momento transferido nulo; sus caractcristicas mas notables son la cxistcncia de funciones dc 
correlacion no analiticas y no univcrsales, y de una singularidad algcbraica en la distribucion 
de momentos sobre la superficie dc Fermi; los modos colectivos bosonicos resultantes poseen 
una relacion de dispersion sin gap [17,40,62], de modo que una cantidad arbitrariamente 
pequeiia de energia es capaz de producir excitaciones a partir del estado fundamental. Sin 
embargo, en situaciones mas realistas se deben tener en cuenta el spin electronico y procesos 
de dispersion electronica mas complejos como dispersion hacia atras y umklapp. En estos 
procesos el momento transferido es 2kp y ikp respectivamente. Estos problemas han sido 
tratados en el contexto de la bosonizacion operational usual en el capitulo 2, siguiendo las 
referencias [17,39]. La inclusion del spin electronico da lugar a la separacion spin-carga, y 
los nuevos terminos de dispersion hacia atras y umklapp generan gaps en los espectros de 
ambos sectores respectivamente, y por supuesto cambian los decaimientos algebraicos de las 
funciones de correlacion por decaimientos exponenciales. 
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En trabajos recientes [35, 66, 67] se trataron estos modelos mas realistas mediante el 
enfoque de la teoria de campos y la bosonizacion funcional. Se introdujo una version no local 
y no covariante del modelo de Thirring en el que las densidades y corrientes fermionicas se 
acoplan mediante potenciales bilocales generales, y una generalizacion a un modelo de Gross- 
Neveu SU(N) no local con N=2 para tener en cuenta el spin electronico. Estos modelos 
contienen al modelo de TL como caso particular (para N=l). Aunque constituyen un marco 
elegante para tratar problemas de muclios cuerpos en ID, presentan dos inconvenientes: i) 
la tecnica de bosonizacion funcional empleada descansa en la regularizacion del jacobiano 
tomada de la teoria de campos covariante, con las consiguientes diferencias halladas con las 
Ref. 17,40,62 mencionadas en el capitulo anterior, y ii) la descripcion de los modelos con 
simetria SU(N) se realizo mediante bosonizacion no abeliana, lo que resulta poco practico 
para el analisis final. 

En el presente capitulo salvamos estos problemas introduciendo un modelo de tipo Thirring 
con dos especies de fermiones, generalizando el modelo de dos fermiones de la Ref. 68, con- 
struido originalmente como una teoria local y covariante, al caso en que las interacciones 
entre corrientes fermionicas se hallan mediadas por funciones bilocales. Este modelo no in- 
cluye interacciones de tipo umklapp. Solo posee interacciones de tipo inversion de spin (IS), 
que en el caso local se reducen a las de dispersion hacia atras. Adicionalmente utilizamos 
la tecnica de bosonizacion abeliana en materia condensada, tal como se la formulo en el 
capitulo anterior. 

El resultado final puede considerarse una extension del trabajo de Grinstein, Minnhagen 
y Rosengren [69] donde se estudio una version simplificada del problema de inversion de spin 
(contenido en nuestro modelo). Extendiendo formulaciones previas de modelos masivos [28, 
32,70], nuestro procedimiento permite obtener una funcional de vacfo bosonizada. El capitulo 
esta organizado de la siguiente manera. En la seccion introducimos el mencionado modelo 
no local de dos fermiones y explicamos su relacion con la descripcion no abeliana previa. En 
la seccion 14.31 establecemos la equivalencia entre la funcion de particion fermionica inicial 
y la que corresponde a una extension no local de dos bosones del modelo seno-Gordon. 
Finalmente, en la seccion 14.41 reunimos los puntos principales de nuestras investigaciones. 

4.2 El modelo y su relacion con una descripcion no 
abeliana previa. 

En esta seccion introducimos una version no local del modelo de Thirring que incorpora el 
spin electronico considerando dos especies de fermiones, donde cada especie representa un 
estado de spin. La accion euclidea inicial es 
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d^xd^y^\x)j,^\x)U(^,){x - y)^\y)j,^\y) (4.1) 
donde a,b = 1, 2, con 1 =|, 2 =|, las corrientes son las corrientes fermionicas usuales 

jI 



■■1 - jI = ^'i,^', (4.2) 

y las matrices V^^^ tienen la forma 

^ Vip.) '^(m) '^(m) + ^{p)J 

donde V^^-, y V"^-^ son los potenciales de acoplamiento relacionados a las constantes g2 y g-i de 
dispersion hacia adelante de Solyom [50]. Estas funciones se vinculan a interacciones del tipo 
corriente-corriente que no invierten el spin. Las funcions f/(^)(x — y) son los acoplamientos 
que gobiernan los procesos que cambian el estado del spin electronico. Hemos dejado las 
constantes g j Qs para facilitar la comparacion con versiones locales del modelo. Por ejemplo, 
el caso fiffi = y V^^^ = V^^^ = 5^(x — y) corresponde a dos modelos de Thirring usuales 
desacoplados. La accion ()4.H1 tiene simetria quiral U(l) manifiesta, es decir, es invariante 
bajo la transformacion ^nsS^a^ _^ ijragnsf conservacion del spin fermionico 

tambien se preserva en esta teoria. 

La teorfa definida mas arriba es similar a la descripta por Zinn-Justin(ZJ) [68]. Existen, 
sin embargo dos diferencias importantes. En primer lugar nuestro modelo tiene en cuenta 
la posible naturaleza de largo alcance de los potenciales, mientras que el modelo de ZJ es 
local. Por otro lado, el modelo de ZJ no incluye terminos del tipo g4 asociados a procesos de 
dispersion fermionica que involucran solo una rama (derecha o izquierda), ni para diagramas 
de IS ni para diagramas ordinarios. 

En lo concerniente a la relacion entre nuestra accion y modelos previos inspirados en 
materia condensada, debemos mencionar los primeros trabajos de Luther y Emery [71], y 
Grinstein, Minnhagen y Rosengren [69]. Los primeros introdujeron el Uamado modelo de 
dispersion hacia atras. Aunque este sistema, en principio, no posee IS, en el limite local 
que ellos consideraron, los diagramas de dispersion hacia atras coinciden con aquellos que 
cambian el spin. Grinstein y col. incluyeron desde el comienzo interacciones de tipo IS 
teniendo en cuenta un potencial coulombiano. Aunque dicho modelo es no local, considera 
el mismo potencial para todas las clases de diagramas (IS y los ordinarios). Por otro lado, 
con el objeto de establecer una relacion entre su teorfa y un sistema de tipo gas de Coulomb, 
los autores consideraron de nuevo el limite local. Y, nuevamente, no incluyeron terminos de 
tipo g^. 
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Mostremos ahora que la acci6n ()4.1|) puede escribirse de un modo alternative. Consider- 
emos las corrientes U(N) 



J'. 



a 



(4.4) 



con 



(4.5) 
(4.6) 



siendo P los generadores de SU(N) normalizados de acuerdo a 



(4.7) 



Con estas corrientes, podemos definir un modelo de Gross-Neveu, invariante quiral y no 
local, con una accion dada por 



S= Sx^^"^- d^xd%J^ix)V^fJx-y)Ji^iy), a, /3 = 0, 1, iV^ - 1 (4.8) 



Este modelo no abeliano fue considerado en la Ref. 35. La accion efectiva bosonizada 
obtenida mediante bosonizacion no abeliana dio lugar a una funcional de Wess-Zumino- 
Witten (WZW), que resulta dificil de tratar a la hora de obtener el espectro fisico. En 
la proxima seccion mostramos como esta tarea se ve simplificada comenzando por ()4.1|1 en 
lugar de ()4.8|) y combinando bosonizacion abeliana en el marco de la integral funcional y la 
aproximacion armonica autoconsistente. 

4.3 La accion bosonica equivalente 

Comenzamos considerando la funcion de particion 





(4.10) 
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donde M es una constante de normalizacion. Es conveniente escribir 



S — Sq-\- (Sflip, 



(4.11) 



donde 



5o = y" d^x^^^'^'^ - ^ 1 <exd^yjl{x)V^;~^{x - y)jj(y) (4.12) 

y 

S^^^-Qsj d'xd^y^\x)^^^\x)U^^){x-y)^\yh^^\y). (4.13) 

La razon para esta separacion se halla en el hecho de que -S'o contiene todos los terminos 
de interaccion que poseen invarianza quiral separada para cada especie fermionica (estados 
de spin). Son interacciones de tipo Thirring, es decir que se transforman en terminos libres 
de la accion. El segundo termino no posee invarianza quiral separada, y sera expandido en 
una serie perturbativa, en analogia con lo hecho con el termino de masa en la bosonizacion 
funcional del modelo de Thirring masivo [32]. 

De igual forma que en el capitulo anterior, la introduccion de campos vectoriales auxiliares 
mediante la transformacion de Hubbard- Stratonovich permite escribir 



- j d^x {0 + g^") - S[A] - 5flip 



(4.14) 



donde H' es una nueva constante de normalizacion que incluye al determinante fermionico 
libre, y 



con 



S[A] = \j d'xd'y {v^-^y\x-y)A;{x)Al{y), 

/ \ab 

f V^~| j definido a traves de la ecuacion 

/ d'y (v^-^y (x - y)V^%{y - z) ^ 5^'\x - z)S^'^. 
Descomponemos ahora yl^ en sus partes transversal y longitudinal 



(4.15) 



A';^{x) = e^,d,r{x) + dX{x), 



(4.16) 



(4.17) 



donde 0" y ry" son campos escalares. Ademas realizamos un cambio en los campos fermionicos 

^"{x) = e-^l^^^^^^^-^^^^^^lx^l^^) (4.18) 



(4.19) 
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cuyo Jacobiano fue calculado en el capitulo anterior. Obtenemos entonces 



siendo Sef una suma de tres partes: 



donde 



'S'ef — 5'oF + •S'oB + 5'flip 



^0F = J d'x (x'^^x'+xW), 



(4.20) 



(4.21) 



(4.22) 



X [t^.t^pd,r{x)dp(t)\y) + d,,7]\x)d^Ti\y) + 2t^,d,r{x)dpv\y)\ , (4.23) 



y S'flip es el mismo termino de interaccion IS ya definido en la Ec. ()4.13|1 . Respecto de este 
ultimo termino, desde ahora nos restringiremos al caso de interacciones IS de contacto: 



U{x- y)(o) = U{x- ?/)(i) = 5^^\x - y), 



(4.24) 



y por medio de una transformacion de Fierz seguida del cambio quiral definido en las Ecs. 
fl4.18|) y ()4.19|) . po demos escribirlo en la forma 



5'fiip = 2gs d X 



(4.25) 

Ahora estamos listos para hacer la expansion de la funcion de particion tomando Qs como 
parametro perturbativo: 



n=0 



' n 

n 



Xi 



ii=l 



2g(</,i -</,2) 1 ]_+ 75 ^ ,1 ^-,2 1 - 75 _2 



-X -x 



. , -,,,i-<a2)_i1 - 75 1 -2! +75 2 
X +e-^- >x — — X 'X ^ X 



OF 

(4.26) 



donde ( )qp significa valor medio en una teoria con accion S'qf- Solo valores medios que 
involucran un mismo mimero de factores de la forma ^'yCiX + 75)%" Y ix"(l ~ 75)x" (^n esta 
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expresion no debe entenderse suma sobre indices repetidos) son no nulos, por lo tanto la 
funcion de particion puede escribirse en la forma 



J n=0 ^'^■^ J \i=l 



><{Ilfi^^)^-T^X^xmyO^-^x'{y^)) . (4.27) 



OF 



El siguiente paso consiste en introducir dos campos escalares locales sin masa 7?" que seran 
asociados a los fermiones libres x° y x". Este truco permite reemplazar los valores medios 
fermionicos en la expresion de arriba por sus contrapartes bosonicos regularizados [32,68], 
lo que conduce a 




(4.28) 



Ahora usamos el hecho de que las divergencias infrarojas del propagador del campo i?" 
proveen una condicion de neutralidad para los valores medios de los operadores de vertice. 
Esto significa que el valor medio (nr=i ^*^''^^^'^)oi? ®^ solamente si — [68]- 

La condicion de neutralidad permite rearmar la scric perturbativa de un modo no trivial. 
Llegamos entonces a la accion completamente bosonizada S'bos 



^bos = ^0B + J (fx l^^{d,rf-^cos[2ig{4>'-<P') + V^{^'+^')\y (4.29) 

En este punto es interesante observar que existe un cambio de variables que permite 
expresar la accion en un modo muy sugestivo. Efectivamente, escribiendo 



(4.30) 
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(4.31) 

(4.32) 
(4.33) 



donde el signo mas (menos) corresponde al par 0^,77^(^^,77^), se ve que el campo 9 se de- 
sacopla completamente de los otros y por lo tanto puede integrarse, y la accion bosonica 
queda 



'S'bos — Sp + Scr, 



donde 



(4.34) 



^p-^J [i^M' - (diVp? - 29i0,9o77,] +\J d'xd'y (v^^))"' {x - y) 

X [e^iye^,pdy(j)p{x)dp(t)p{y) + dp,r]p{x)dp,r]p{y) + 2€p,^d^(f)p{x)dp,r]p{y)] , (4.35) 



S-^^J d'x [(9i0,)^ - (9i77,)^ - 2910,9077,] +^ J d'xd'y {V^^^^y' (x - y) 
X [enuep.pd„(pa{x)dp(j)aiy) + dp,7]„{x)dp,r]a{y) + '2e^^d^(l)„{x)dp,r]^{y)] 



I 



+ / d^x 



con las funciones ( V^^^* 



definidas como 



(4.36) 



(4.37) 



Esto, a su vez, conduce a una factorizacion de la funcion de particion en la forma Z — 
ZpZfj. Este resultado es una clara manifestacion de la separacion spin-carga [69,71]. Zp 
es la funcion de particion asociada a las excitaciones de densidad de carga. Coincide con 
el modclo dc Tomonaga-Luttinger sin spin estudiado en el capitulo anterior y en la Ref. 
41. Zfj describe excitaciones de densidad de spin. La accion 5", corresponde a un modelo 
seno-Gordon con termino cinetico no local, similar al considerado previamente en la Ref. 70. 
En el capitulo siguiente derivaremos una expresion para el gap de su espectro en funcion de 
los potenciales V^^^. 
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Dado que nuestro objetivo principal es analizar el sector de spin, desde ahora enfocaremos 
nuestra atencion en S^- Tranformando Fourier la accion (con la excepcion del termino del 
coseno, cuya transformada de Fourier no es muy esclarecedora) se obtiene 



(p{p)(p{-p)A{p) + fj{p)fj{-p)B{p) + <i){p)f}{-p)C{p) + —9{p)e{-p) 



con 



- j d?x cos [y^ig(t)a + VStt^) (4.38) 



C =piPo 



(4.39) 
(4.40) 
(4.41) 



En la expresion dc arriba, <j){p), T){p) j 6{p) son las transformadas de Fourier de (j)a{x), 
V(t{x) y 9{x) respectivamente. Es conveniente considerar todavia otro cambio de variables 
que diagonaliza la parte cuadratica de Sa- Este cambio esta dado por 



e 



An 



A7r + 25c/V \/7r 



i^gCp' 



(4.42) 
(4.43) 

(4.44) 



Att + 25^2^2 

donde C Y X son los nuevos campos bosonicos, y A = — 4AB. La accion resultante se 
lee 



+ 



(fp 

(27r)2 2K, 
d?p 



— (Po + ^.P?)e(p)e(-p) fd'x cos(V8^0 

C(P)C(-P) 



2. p4(^ _ ^2^. 



+ 



(27r)2 2{plnV([^+plV^o){7i-g^V^i) 
(fp 1 



(27r)2 2 



— + 



(0) 



TT 



^{p)^{-p) 



(4.45) 
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donde 




1 _L r_ya 



^ ^(1) 



(4.46) 



(4.47) 



Se observa que los campos C y se desacoplaron completamente de ^. La parte de la 
accion que depende de ^ coincide con la obtenida mediante bosonizacion operacional en la 
seccion 14.31 y es la parte relevante a considerar. Mas generalmente la reescribimos como 



Scr 



con 



(27r)2 
F{p) 



j d^x cos(V8^0, (4-4J 



(4.49) 



Esta accion corresponde a un modelo seno-Gordon con termino cinetico no local, introducido 
en la Ref. 70. 



4.4 Conclusiones 

En este capftulo hemos mejorado una version no local del modelo de Thirring que provee 
una descripcion tratable de los Ifquidos de Luttinger con spin, basada en teon'a de campos. 
Efectivamente, en el contexto de teorias del tipo Thirring no local, los tratamientos previos 
de las interaccciones de inversion de spin condujeron a un complicado modelo no abeliano 
(ver por ejemplo [35]). Especificamente construimos una accion basada en dos especies de 
fermiones que permite tener en cuenta interacciones de inversion de spin de un modo elegante 
y simple. Aunque nuestro modelo se inspire en el considerado en la Ref. [68], incluye inter- 
acciones no contenidas en ese trabajo previo (los llamados diagramas en la terminologia 
de Solyom [50]). Ademas, la teon'a que presentamos posee potenciales bilocales generales 
que gobiernan las interacciones que no invierten el spin. Parametrizamos estos potenciales 
en terminos de las funciones y que se asocian a la dinamica de las densidades de 
carga y spin respectivamente, una vez que la separacion spin-carga se liace manifiesta luego 
de un cambio de variables apropiado (ver ecuacion ()4.32|1 y ()4.33j) ). Aunque nuestro analisis 
es unicamente valido para interacciones de inversion de spin locales, pudimos mantener la 
dependencia con la distancia en los potenciales ordinarios (que no invierten el spin) has- 
ta el final de los calculos. Bajo estas condiciones obtuvimos una accion bosonica efectiva 
cuyos grados de libertad de carga coinciden con descripciones previamente encontradas en 
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el capitulo anterior del problema de dispersion hacia adelante sin spin. Respecto del sector 
de spin, que es el de mayor interes en el presente contexto, hallamos que se corresponde con 
un modelo seno-Gordon con termino cinetico no local. En el capitulo siguiente estudiaremos 
su espectro dentro del marco de la aproximacion armonica autoconsistente, derivando una 
expresion para el gap de las excitaciones de baja energia. 
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Capitulo 5 

La aproximacion armonica 
autoconsistente 



Examinamos la aproximacion armonica autoconsistente y su iitilizacion en el 
marco de la integral funcional y las tcorias clc materia conclcnsada. En particular 
hallamos una formula para el gap de las excitaciones del sector de spin del 
modelo estudiado en el capitulo anterior, como funcion de potenciales arbitrarios 
de interaccion electron-electron. Proponemos ademas un nuevo metodo para 
determinar cl parametro incognita asociado a esta aproximacion. Comprobamos 
la validcz de esta nueva tecnica en el contexto del modelo seno-Gordon y como 
aplicacion no trivial consideramos el regimen de escala del modelo de Ising en 
2D fuera del punto critico y en presencia de un campo magnetico h. En este caso 
derivamos una expresion aproximada que relaciona la longitud de correlacion ^, 
T — Tc y h. Los resultados de este capitulo constituyen aportes originales de 
esta tesis [42, 43]. 



5.1 Introduccion 

La aproximacion armonica autoconsistente (SCHA, self-consistent harmonic approximation) 

es una tecnica no perturbativa que ha sido utihzado extensivamcntc en aphcaciones de 
mecanica estadistica [72,73] y materia condensada [43,74-77]. Consistc en reemplazar una 
accion vcrdadera S por una accion de prucba 5*0 que hace que cl problcma resulte tratable. 
Usualmente 5*0 es una accion cuadratica que depende de cierto parametro desconocido fl. 
Este parametro debe ser determinado mediante algun criterio como por ejemplo la mini- 
mizacion de la energfa libre. Esta aproximacion esta mtimamente relacionada al potencial 
efectivo gausiano [78-80] en teoria cuantica de campos, una aproximacion variacional al 
potencial efectivo que utiliza como cstado fundamental de prucba una funcional dc onda 
gaussiana dependiente de un parametro de masa. Tambien esta basado en un principio de 
minima sensibilidad [81,82] para determinar el parametro adicional. 
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En el capitulo anterior mostramos que la inclusion de terminos que invierten el spin en 
la accion de una teoria de electrones fuertemente correlacionados en una dimension espacial, 
condujo a una accion bosonica de tipo seno-Gordon para las fluctuaciones de los grados de 
libcrtad de spin. Aunque el espectro del modelo seno-Gordon usual se conoce exactamente 
a partir de los trabajos de Dashen, Hasslacher y Nevcu (DHN) [83], en el caso que nos 
ocupa la presencia de interacciones de largo alcance dan lugar a un termino cinetico no local 
que destruye la solubilidad. Esta situacion nos Ueva naturalmente a considerar metodos 
aproximados. En el presente capitulo nos proponemos hallar el espectro de excitaciones 
del mencionado sector mediante la SGHA. El prcscntc calculo resulta de interes por dos 
motivos adicionales: en primer lugar, la SCHA es mucho mas sencilla de implementar que 
las tecnicas semiclasicas utilizadas en los mencionados trabajos. Y en segundo lugar, la 
SGHA puede extenderse facilmente a modelos mas complicados como por ejemplo el doble 
seno-Gordon [84] en teon'a de campos, o la version continua del modelo de Hubbard extendido 
en Uenado medio [85] y Ifquidos de Luttinger acoplados [86]. 

Si bien nuestra motivacion inicial para abordar el estudio de la SCHA fue contar con 
un metodo de aproximacion para poder estimar el valor del gap (en modelos en los que la 
no localidad impide hallar una solucion exacta) al familiarizarnos con la tecnica pudimos 
realizar algunos aportes originales en lo concerniente al metodo en si mismo. Esto nos 
ha Uevado a presentar tanbien en este capitulo, a modo de leve digrcsion, los detalles de 
nuestra propuesta. En particular, senalamos que en problemas en dos dimensiones existe 
un modo alternativo de determinar el parametro fl. Este metodo esta basado en teoria de 
campos conformes [87,88]. Mas aiin, mostramos que nuestro metodo conduce a mejoras en 
los resultados para el modelo SG respecto del SGHA estandar y nos permite dar una nueva 
descripcion del modelo de Ising bidimensional (MI2D) fuera del punto cn'tico. En el primer 
caso explotamos la existencia de resultados exactos [83,89] para verificar la consistencia de 
nuestra propuesta, obteniendo una respuesta cualitativamente buena para la masa del soliton. 
Aplicamos entonces las mismas ideas al MI2D a T 7^ Tc y 7^ 0, un modelo no integrable 
para el cual se conocen pocos resultados cuantitativos [90-92] . Utilizamos la representacion 
fermionica del MI2D. Dado que la SGHA estandar esta restringida a modelos bosonicos, el 
nuevo procedimiento provee tambien una extension de la aproximacion gausiana a teorfas 
fermionicas bidimensionales. Nuestro resultado principal es una ecuacion algebraica que 
permite obtener el comportamiento de la longitud de correlacion como funcion de T — Tc y 
h. 

Dcbcmos cnfatizar que no estamos introducicndo una nueva aproximacion, sino solo un 
metodo para determinar el parametro. Como es bien sabido, la SCHA es una aproximacion 
no controlada, es decir, no hay ningiin parametro pcrturbativo involucrado. Es claro entonces 
que la misma critica puede hacerse a la presente propuesta. 
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5.2 Detalles de la aproximacion y el espectro del sector 
de spin en electrones unidimensionales 

Comenzaremos mostrando el desarrollo mediante integrales funcionales de la SCHA. En 
general comenzamos con una funcion de particion 

Z= [ Vfie-^ (5.1) 



donde Vfi es una medida de integracion funcional bosonica. Una manipulacion elemental 
conduce a 

Z = J ^j^^^_s. j Vf^ e-^o = 2, (e-(^-^")>„ . (5.2) 
para cualquier accion de prueba Sq. Por medio de la propiedad 

(e-^) > e-^f\ (5.3) 

para / real, y tomando logaritmo natural en la ecuacion ()5.2p . obtenemos la desigualdad de 
Feynman [93]. 

\nZ>\nZo- {S - So)o (5.4) 

La aproximacion consiste en reemplazar la accion verdadera, difi'cil de tratar, por una 
accion de prueba mas simple, que contenga algun parametro libre. Este se fija maximizando 
el lado dereclio de la desigualdad ()5.4p . 

Consideramos a continuacion como accion verdadera, la accion que describe la dinamica 
del sector de spin en teorfas de electrones altamente correlacionados en una dimension ()4.48|1 . 
obtenida en el capitulo anterior. Ademas como accion de prueba tomamos una accion 
cuadratica de la forma 



2 



(27r)2 

Para la situacion mencionada, F se define como 



^(p)^^(-p) + Ye(p)e(-p) 



(5.5) 



F{p) = -^{pI + v.pI), (5.6) 

J^sVs 

aunque para los desarroUos siguientes consideraremos una F arbitraria, y obtendremos re- 
sultados generales. Una vez realizada la sustitucion de acciones es inmediato obtener el 
espectro: 

F{p) + = 0. (5.7) 
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Volviendo a frecuencias reales, po = iu, Pi = k, j tomando el termino cinetico ()5.6p se 
obtiene la siguiente ecuacion 



VsKsVl'^ + vlk"^ - = 0. 



(5i 



Como se dijo anteriormente, el parametro VL puede determinarse maximizando el lado 
derecho de la ecuacion ()5.4j) . Para alcanzar este objetivo primero escribimos 



In Zq = \q. I T>^ exp 



cPxi{x)(Ai){x) 



- tr InA + const, 
2 



[ln(deti)"i/2] + const (5.9) 

(5.10) 

donde el operador A esta definido, en espacio de Fourier, por 

{Ai){p) = [F{p) + n']i{p). (5.11) 



Es facil entonces obtener 



trlnA = V 



(fp 
(2^ 



\n[F{p) + 



(5.12) 



donde V es el volumen (infinito) de todo el espacio / (fx. Por otro lado, {S — Sq) se puede 
calcular directamente siguiendo, por ejemplo, los pasos explicados en la Ref. 70 y la identidad 
del apendice O El resultado es 



{S - So)q = d X 



TC 



V — — exp 



-An 



(Pp 1 
(27r)2 F(p) + ^2 



+ V- 



f (Pp 



2 J (27r)2 F{p) + D?' 



(5.13) 
(5.14) 



Finalmente, podemos reunir todos los terminos y escribir 



\nZ^-{S-S^)^ = V 



^exp (-47r/o(^])) + ^/o(^^) - \h{n) 



+ const. (5.15) 



donde definimos las integrales 



(2^ 
(fk 



ln[F(p) + 
1 



(27^)2[F(p)+^]2]n+l 
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(5.16) 
(5.17) 



con las propiedades formales 



'-T-'^'^ (5.18) 

^^^ = -2{n + l)ni.^.,. (5.19) 

Hallando el extreme de la expresion ()5.15j) con respecto a f2, y asumiendo que es 
no nulo (una condicion que vale para la mayoria de los potenciales realistas), finalmente 
obtenemos la ecuacion del gap. 

^2 _ ^^-4./„(n) ^ Q_ (5 20) 

La ecuacion ()5.2Up es uno de los resultados principales de este capitulo. Dentro de la 
aproximacion armonica autoconsistente, da una expresion cerrada para el gap como funcional 
de los potenciales V(^-)(p). Resulta interesante obtener un valor explicito para Q para el 
potencial de contacto, dado por 

^(o)(p) = V^d)(p) = l- (5.21) 

En este caso Io{^l) es infinita. Si utilizamos la misma regularizacion empleada en el capitulo 
anterior se obtiene 



A2 7r2 



1 + 



(5.22) 



Esta ecuacion puede facilmente resolverse para A ^ i7 y i7 ^ A. Los resultados estan dados 
respectivamente por 



,2\ -9^ A 
TT / \ TT 



+ ^ (5.23) 



= (5.24) 



5.2.1 El modelo seno-Gordon en teoria de campos 

Si bien el modelo estudiado en la seccion precedente coincide con el modelo seno-Gordon 
considerado tradicionalmente en la literatura de la teoria de campos en el caso local dado 
por ()5.2H1 . existe una diferencia en el tratamiento de las divergencias en las integrales. 
Mientras que en teoria de campos se busca implementar mecanismos que hagan finitos los 
resultados finales, lo que se consigue al renormalizar los parametros de la teoria, en materia 
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condensada las divergencias que aparecen son artificiales, es decir son introducidas por las 
tecnicas mediante las cuales se trata la teoria. Por ejemplo, en el modelo de Hubbard original, 
el espectro es acotado, es decir, las sumas o integrales sobre k se realizan en una region finita. 
Al considerar la linealizacion de la relacion de dispersion, y la inclusion de infinitos estados 
con energfa por debajo del nivel de Fermi, el modelo se vuelve tratable matematicamente 
mediante bosonizacion, pero al costo de la introduccion de divergencias ultravioletas. Por 
otro lado, el parametro que se introduce para regular el modelo [la constante a en la ecuacion 
()2.22|1 . o A en ()4.28j) ]. usualmente tiene un significado fisico: esta asociada con la constante de 
red. Por esta razon es usual que se deje su dependencia explfcita en funciones de correlacion 
y otras cantidades de interes. 

En la seccion anterior seguimos este ultimo camino y obtuvimos una expresion para el gap 
que dependia expllcitamente del cutoff A. En esta seccion seguiremos el camino usual de la 
teoria de campos, es decir, introduciremos una constante de acoplamiento renormalizada, y 
obtendremos expresiones para el gap del modelo seno-Gordon que dependen de ella. Partimos 
de la accion 

donde f{p) es un campo escalar y a es la constante de acoplamiento renormalizada^ mediante 
un orden normal [28,32]; p es el parametro que implementa el orden normal. Por simplicidad, 
en esta formula hemos escrito el termino cinetico en espacio de Fourier, pero hemos dejado 
el termino de interaccion en espacio de coordenadas. 

Como accion de prueba, proponemos una accion cuadratica, 



5*0 



2 



F(p) Q 



(5.26) 



(2vr) 

donde Q es el paramtro de prueba. Con el objeto de realizar el proceso de minimizacion 
estandar, primero evaluamos (5* — 5*0). El resultado es 



{S-So), = V 



2 



p z 



(5.27) 



Insertando ahora ()5.27|) en la desigualdad de Feynman ()5.4|) . y minimizando el lado 
derecho de dicha ecuacion con respecto a Q, finalmente obtenemos 

Esta ecuacion del gap permite extraer una respuesta finita para Q, dependiente del parametro 
de masa p (la diferencia /o(^^) — hip) es finita). Notese que el valor de p es completamente 

^Aqui el termino renormalizado es un abuso de lenguaje, ya que en modelos en (1+1) dimensiones con 
interacciones sin derivadas, el orden normal elimina completamente todas las divergencias 
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arbitrario, si uno lo elige igual a la masa de prueba f2, la solucion de la ecuacion es 



fi^ = a. (5.29) 
El mismo resultado se obtiene si en lugar de p = se toma p = y/a. 

5.3 Determinacion de Q mediant e tecnicas de teoria de 
campos conformes 

Present aremos a continuacion una ruta alternativa para determinar Q. Para este fin ex- 
plotaremos una prediccion cuantitativa de la invarianza conforme para sistemas en 2D 
en el regimen de escala, fuera del punto cn'tico. Partiendo del Uamado 'Teorema-c' [94], 
Cardy [95, 96] mostro que el valor de la anomalia conforme c, que caracteriza al modelo en 
el punto critico, y el segundo momento del correlador de densidad de energia en el regimen 
de escala de la teoria no crftica estan relacionados por 

d'x \x\' {e{x)em = 3 ^2 - A,r ^^'^^^ 

donde e es el operador de densidad de energia, es su dimension de escala y t oc (T — Tc) es 
la constante de acoplamiento del termino de interaccion que saca al sistema fuera del punto 
critico. La validez de esta formula ha sido verificada explfcitamente para varios modelos 
[95,96]. Para el modelo SG, el operador de densidad de energia esta dado por el termino del 
coseno, su dimension conforme es = /5^/47r, t es la constante de acoplamiento a/P"^ y la 
teoria conforme bosonica fibre asociada posee c = 1. 

Nosotros afirmamos que fl puede ser determinado de una forma no variacional comple- 
tamente diferente, forzando la validez de la identidad conforme de arriba para la accion de 
prueba. En otras palabras, impondremos que se verifique la siguiente ecuacion: 

a'^ f 1 

—7 / (i^xIxP (cos 3ip(x) cos 6(p(0))r, = — , (5.31) 

(34 J I I \ /^^v ; /-^v j/o 3^ (2-g)2' ^ ' 

que debe ser considerada como una ecuacion para el parametro de masa Q. Por supuesto, si 
uno esta interesado en comparar la respuesta dada por esta formula con el resultado dado 
por la SCHA usual, al evaluar el lado dereclio de ()5.31|) . se debe adoptar una prescripcion de 
regularizacion equivalente al orden normal implementado en la SCHA. Un calculo cuidadoso 
conduce a la siguiente ecuacion del gap: 

i-r-' = l,^-^ (5.32) 
p p^ 62 

donde hemos definido la variable u = jS'^/An (0 < u < 2) y p es el parametro de orden normal, 
como antes. Vemos que, igual que en la ecuacion SCHA estandar ()5.28|1 . hay diferentes 
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respuestas para diferentes elecciones de p, pero en este caso, los resultados obtenidos para los 
valores ^Ja y f2 son diferentes. En cualquier caso se obtiene una dependencia no trivial de f2 
con en contraste con la SCHA variacional. Esto es interesante si recordamos el significado 
fisico de la masa del gap en el contexto del modelo SG. Efectivamente, como es bien sabido, 
Daslien, Hasslacher y Neveu (DHN) [83] calcularon mediante tecnicas semiclasicas el espectro 
de masa para el modelo SG. Consiste en un soliton (asociado al fermion del modelo de 
Thirring) con masa 



con N = 1,2,... < (2 — u)/u. (De esta ultima condicion es facil ver que para tener 
estados ligados debemos tener u < 2/(A^ + 1). Como consecuencia no hay estado ligado 
para u > 1). Mas recientemente, Zamolodchikov [89], reinterpretando resultados obtenidos 
mediante el ansatz de Bethe, dio expresiones exactas para este espectro. En particular para 
el soliton su formula acuerda muy bien con ()5.33|) . excepto para u cercano a 2, donde predice 
una divergencia. Por simplicidad, aqui comparamos nuestros resultados con la ecuacion 
()5.33|1 . Lo primero que debe notarse es que las masas en el espectro del modelo SG tambien 
dependen de u, igual que nuestra prediccion dada por la ecuacion ()5.32|) . Por lo tanto, al 
respecto, nuestra propuesta parece mejorar la prediccion gausiana estandar para el modelo 
SG, al menos cualitativamente. Para efectuar una discusion cuantitativa mas especi'fica 
comparemos las ecuaciones ()5.32|1 y ()5.33j) como funciones de u. Fijamos p = -y/a, lo que 
corresponde a la prescripcion empleada por DHN al derivar ()5.33j) y ()5.34j) . El resultado 
se muestra en la figura 15.11 donde se puede observar una analogia cualitativa general entre 
ambas curvas. En particular, para 0.7<n<l(n = l corresponde al punto del fermion libre 
del modelo de Thirring y al punto de Emery en el modelo de dispersion hacia atras [71]) 
nuestra prediccion esta de acuerdo con los valores de la masa del soliton calculados por DHN. 
Queremos recalcar que para u = 1 obtenemos ^/^/a = ■\/3/32 ^ 0.30 mientras que el valor 
dado por ()5.33j) es l/vr ^ 0.31 (la SCHA estandar da, por supuesto, fl/y/a = 1). 

5.3.1 El modelo de Ising en 2D 

Habiendo verificado la admisibilidad de nuestra propuesta en un modelo en el que se conocen 
resultados exactos, es ahora deseable explorar un problema no trivial. Consideremos el 
modelo de Ising en 2D fuera del punto critico {T ^ Tc j h ^ 0): 




2-u 



(5.33) 



y una secuencia de estados ligados con masas 




(5.34) 




(5.35) 
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1 - 



0.8 - 



0.6- 



0.4 - 



0.2 - 




0.5 1 1.5 2 



Figura 5.1: Masas en unidades de ^/a como funcion de u. La Imea punteada es Mgoi/v^? 
mientras que la Imea Uena representa ^j^fa. 



donde S'm es la accion critica, t oc {T — Tc), y e{x) y a{x) son los operadores de densidad 
de energi'a y spin respectivamente. Usaremos la representacion fermionica de la mencionada 
accion. De este modo, S'm es la accion de un fermion de Majorana libre y sin masa, y 
e oc Por otro lado, la expresion de (t{x) en terminos de los campos de Majorana es 

mas complicada. Efectivamente, por medio de una transformacion de Jordan- Wigner puede 
escribirse como exponencial de un bilineal fermionico. En analogfa con la SCHA usual, 
proponemos la siguiente accion cuadratica de prueba 

5'o = Sm ~^ ^ J '^^3;e(x), (5.36) 
La ecuacion conforme ()5.3n|l para el presente caso toma la forma 
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ciV [t^ (2 







+ 2t/i(2- A,) (2- A^) (6(r)a(0))J 



1 



(5.37) 



6 TT 



donde fijamos c = 1/2, que es la carga central de los fermiones libres de Majorana, y 
Ag = 1 y Ao- = 1/8 son las dimensiones de escala de los correspondientes operadores. 
Ahora debemos evaluar el valor medio en la accion de prueba. Esto nos dara una ecuacion 
para f2 como funcion de t y /i. Las funciones de correlacion energia-energia y energia-spin 
fueron calculadas por Hecht [97] mientras que el correlador spin-spin puede hallarse en el 
trabajo de Wu, McCoy, Tracy and Barouch [98]. Como es usual, se define la longitud de 
correlacion ^ = l/4f2 y se considera el limite de escala, dado por ^ — > cxd, r — > cxd, con r/^ 
fijo. El siguiente paso es usar la expresion de los correladores para {r/S,) « 1 y realizar 
las integrales correspondientes. En este punto debemos tener en cuenta que las funciones de 
correlacion son proporcionales a ciertas funciones de escala F±{r/^) donde los signos + y — 
corresponden a los casos f2 > y f2 < respectivamente. En otras palabras, el parametro 
fl puede verse como un parametro que define una nueva temperatura cn'tica efectiva. Las 
funciones F± describen el regimen de escala por encima y por debajo de esta temperatura. 
Dado que estamos aproximando una perturbacion magnetica en el sistema, es claro que 
debemos usar la funcion F_. De este modo obtenemos la siguiente ecuacion que relaciona ^, 
hyt: 



donde hemos introducido las constantes numericas Ci = 0.749661 y C2 = 0.186966. El valor 
absoluto del campo magnetico en el segundo termino proviene del heclio de que (ecr) oc {a) y 
el producto (a) h debe ser positivo dado que la magnetizacion y el campo magnetico tienen 
que tener la misma orientacion. Para ^ fijo, esta ecuacion da una dependencia simple de h 
como funcion de t. Efectivamente, para h > tenemos una semi elipse levemente rotada en 
el piano h — t superior, y para h < tenemos su refiexion sobre el eje t = 0. 

Si reescribimos la ecuacion 15.381 en terminos de la longitud de correlacion a campo nulo, 

__8_ 

.^0 = l/4t, y la combinacion adimensional x = I ^ I /4^o obtenemos 



Los signos + y — en el tercer termino del miembro izquierdo corresponden al caso t > y 
t < respectivamente. La accion ()5.35|) define una familia de teorias de campos dependientes 
de un parametro x [91]- 

Con el objeto de verificar la consistencia de la ecuacion de arriba, consideramos los 
Kmites h—^Ojt-^Oen forma separada. El primer caso corresponde a x ^ 00 y se obiene 



t\40^ + Cih\AO^^/^ + C2t|;i|(40'^/^ = 1 



(5.38) 




(5.39) 
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Figura 5.2: Longitud de correlacion en unidades de como funcion de x, para ambos casos, 
t > y t < 0. 



inmediatamente ^ = como era de esperar. En el segundo caso tenemos x — > y obtenemos 

I- — 

entonces ^ ~ | /i | , lo cual esta en acuerdo con el resultado exacto obtenido en las Ref.99, 
100. Debemos mencionar que en estas referencias la constante de proporcionalidad fue 
determinada exactamente en el valor 4, 4, mientras que nuestro calculo aproximado conduce 
al valor 3, 7. Volviendo al caso general, resolvimos la ecuacion ()5.38j) numericamente para 
^ como funcion de x para t > y t < 0. Los resultados se muestran en el grafico de la 
figura 15.21 En el caso t > la longitud de correlacion se incrementa en forma monotona 
desde cero y alcanza el valor correspondiente a. h = 0, desde abajo cuando x — ^ cxd. En el 
caso t < 0, aunque el comportamiento de ^ parece muy similar al caso previo, una mirada 
cuidadosa muestra que presenta una diferencia sutil, mostrada en la figura Para x ~ 2 
la longitud de correlacion pasa sobre el valor ^o, alcanza un maximo y entonces tiende a 
^0 desde arriba cuando x oo. Como este comportamiento depende de los valores de las 
constantes Ci y C2 no sabemos si es efectivamente una propiedad del modelo de Ising o un 
artificio introducido por nuestra aproximacion. 
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Figura 5.3: Igual que la Fig. 15.21 Al agrandar la escala se observan los detalles del compor- 
tamiento de ^(x) para t > y t < 0. 



5.4 Conclusiones 

En este capitulo hemos reconsiderado el conocido metodo de aproximacion armonica auto- 
consistente, en el cual una accion comparativamente compleja es reemplazada por un sistema 
cuadratico mas simple dependiente de un parametro de masa fl que se determina usualmente 
mediante un calculo variacional. Aplicamos este metodo para la obtencion del gap del espec- 
tro de excitaciones del modelo de Thirring no local con interacciones de inversion de spin, 
formulado en el capitulo anterior. 

Al trabajar con la SCHA advertimos que, para el caso de teorfas 1+1 dimensionales, el 
parametro fl se podia determinar de un modo alternativo, no variacional. Nuestra propuesta 
se basa en una consecuencia del teorema-c de Zamolodchikov [94] derivada por primera vez 
por Cardy [95]. Ilustramos la idea considerando el modelo seno-Gordon. Mostramos que 
para este modelo nuestro metodo da una prediccion bastante buena para el comportamiento 
de la masa del soliton como funcion de (ver las ecuaciones ()5.32p y ()5.33p y la Fig. 15. 1|) . 

Como aplicacion no trivial consideramos el modelo de Ising en 2D fuera del punto cn'tico 
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(T 7^ Tc and h ^ 0). A partir de una descripcion mediante teoria de campos en terminos 
de fermiones de Majorana, proponemos una accion cuadratica de prueba dependiente de un 
parametro Q que define una longitud de correlacion aproximada ^. Nuestro resultado prin- 
cipal esta dado por la ecuacion ()5.H8j) (o su forma alternativa ()5.39j) ) que permite determinar 
el parametro Q (o lo que es lo mismo, ^) en terminos de los parametros ffsicos originales t y 
h. 

Seria interesante probar nuestro enfoque en otros modelos tales como la version continua 
del modelo de Ising tricritico, que puede describirse mediante el segundo modelo de la serie 
minimal unitaria [101,102] con carga central c = 7/10. 
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Capitulo 6 

Interacciones spin-orbita 



En este capitulo calculamos funciones de corrclacion en sistenias unidiniension- 
ales de electrones en interaccion en los que los grados de libertad de carga y 
spin se encuentran acoplados a traves de la Interaccion spin-orbita. Estudiamos 
fluctuaciones de tipo ondas de densidad de carga y spin, y de tipo supercon- 
ductor singulete y triplete. Mostramos que la interaccion spin-orbita niodifica 
los exponentes del decainiiento de las funciones de corrclacion y el diagrama de 
fases del sistema. Ademas encontramos que susceptibilidades que eran finitas a 
bajas temperaturas, se vuelven divergentes a causa de la interaccion spin-orbita. 
Estos resultados constituyen una contribucion original de esta tesis [44] 

6.1 Introduccion 

Al considerar el comportamiento de los electrones dentro de materiales, debe tenerse en 
cuenta que estos se mueven en presencia de campos electricos. Como consecuencia, ex- 
perimentan no solo la fuerza electrostatica originada en estos campos, sino una influencia 
relativista conocida como interaccion spin-orbita (SO) que rompe la simetria de rotacion de 
spin SU(2). Su origcn sc cncuentra en el acoplamicnto dc Pauli cntrc cl momento magnetico 
de spin del electron y un campo magnetico que aparece en su sistema de referenda en reposo 
debido al movimiento en un campo electrico. Una forma general de describir la interaccion 
SO consiste en agregar el siguiente termino al Hamiltoniano, que se obtiene a partir de la 
expansion cuadratica en v/c de la ecuacion de Dirac: 



Aqui mo es la masa en reposo del electron, p es el operador momento, a — {a^, Cy, a^} es el 
vector dc las matrices de Pauli, y V{r) es el potencial de la particula. 

En materiales cristalinos tridimensionales, la energi'a V^(r) proviene exclusivamente del 
potencial cristalino micro scopico. Dresselhaus [103] mostro que en estructuras cristalinas 
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sin simetria de inversion como la estructura de tipo zinc-blenda (por ejemplo el material 
semiconductor GaAs posee esta estructura), la interaccion SO conduce a un splitting de la 
banda de conduccion en dos subbandas. La magnitud del splitting es proporcional al cubo 
del numero de onda k del electron. 

Es posible obtener gases de electrones bidimensionales partiendo de sistemas en 3D con- 
finando el movimiento electronico a dos dimensiones mediante la aplicacion de un campo 
electrico perpendicular, generando un pozo cuantico. La reduccion de la dimension efecti- 
va disminuye la simetria del cristal subyacente, y agrega un termino adicional, lineal en k, 
al splitting de las subbandas. Mas aiin, si el pozo cuantico es suficientemente estrecho, el 
termino lineal en k se vuelve dominante. Por otro lado Rashba [104] mostro que este potencial 
confinante macroscopico da lugar a un segundo termino en el Hamiltoniano de interaccion SO 
que es responsable de la aparicion de un termino lineal en k adicional en el splitting de las sub- 
bandas. En diversos sistemas, como lieteroestructuras de InGaAs/InAlAs [105], GaAs [106] 
y GaAs/AlGaAs esta se convierte en la contribucion mas importante al acoplamiento SO. 
En este ultimo caso, el splitting de las subbandas tambien depende del potencial confinante, 
hecho que ha sido comprobado experimentalmente por ejemplo en las Ref. 105-108 mediante 
el estudio de fenomenos de transports 

Las investigaciones tendientes a la comprension y control de los fenomenos que involucran 
el spin electronico en materiales semiconductores ban cobrado un renovado interes en los 
liltimos alios a partir de la idea de fabricar unidades spintronicas, en las que se utiliza el spin 
electronico en lugar de su carga para el manejo y almacenamiento de informacion. En ellas, 
la interaccion SO, y el efecto Rashba en particular juegan un rol central [109-112]. 

Existen diversas tecnicas para crear sistemas de electrones en una dimension a partir 
de gases de electrones bidimensionales [113]. En esencia, todos utilizan un confinamiento 
adicional mediante un potencial transversal. Este se convierte asf en una fuente adicional 
de potencial macroscoopico, origen de la interaccion spin-orbita. Si el confinamiento es 
suficientemente fuerte (angosto y profundo), entonces este campo se vuelve importante frente 
al efecto Rashba, y puede convertirse en dominante. Aunque hasta nuestro conocimiento no 
existe evidencia experimental de la interaccion SO que resulta de este potencial, estudios 
teoricos indican que afecta cualitativamente el comportamiento de la energia de splitting 
como funcion del vector de onda k. En efecto, Moroz y col. resolvieron la ecuacion de 
Shrodinger que resulta de situar un electron estrictamente bidimensional (en el piano xy) en 
un campo electrico transversal (en la direccion del eje z) y en un potencial cuadratico en el eje 
X [114]. Ademas consideraron los terminos de interaccion SO del tipo de la ecuacion ()6.1|) con 
los potenciales anteriormente descriptos. Dado que el sistema continua teniendo invarianza 
traslacional en el eje y, ky sigue siendo un buen numero cuantico. Las energi'as como funcion 
de ky que se encuentran al resolver la ecuacion de Schrodinger se muestran en la figura 16.11 
Se observa un splitting entre las subbandas de spin para arriba y spin para abajo, y ademas 
un nueva caracten'stica, propia del confinamiento a una dimension, que es la deformacion 
de cada una de las bandas como funcion de ky. La caracten'stica mas importante de esta 
deformacion consiste en que cada banda pierde su eje de simetria vertical, y la velocidad de 
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Figura 6.1: Espectro de energia de los electrones en un alambre cuantico con interacciones 
spin-orbita. 



Fermi de los electrones se vuelve diferente para cada direccion de movimiento. Los calculos 
indican [114] que la diferencia de las velocidades de Fermi se incrementa en forma monotona 
con el acoplamiento SO Uegando a ser del orden de 10-20%. Esta razon se encuentra dentro 
del rango en el es que es posible realizar mediciones experiment ales. 

6.2 Formulacion del modelo y formalismo 
6.2.1 Hamiltoniano fermionico 

Como fue mencionado ya varias veces a lo largo de esta tesis, un ingrediente central de los 

electrones en una dimension cn mctalcs y materiales semiconductores es la alta correlacion. 
Dc modo que se vuelve ncccsario incorporar la interaccion elcctronica en modelos para in- 
teracciones SO. Para ello fue propuesto el siguiente Hamiltoniano [37,38] 



H — 'Ho + Hint, 



(6.2) 



donde el Hamiltoniano libre es 
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ipl^ crea un fermion de spin a que se mueve a la derecha (r = +1) o a la izquierda (r = — 1). 
Este Hamiltoniano es similar al correspondiente al modelo de Tomonaga-Luttinger con spin, 
con la diferencia de que en este caso, los fermiones que se mueven a izquierda y derecha en 
ambas bandas, poseen diferentes velocidades de Fermi {vi ^ V2), reflejando la asimetria del 
espectro libre. El termino de interaccion describe interacciones de dispersion hacia adelante, 
y posee la forma usual: 

Hint = dx {g2\\SaP + g2±Sa,-f3) 'ip[,a^R,a1pR^,3lpL,(5 

a, 13 

+ j dx {gi\\5c,p + 5'4±5„ V'i„^r,a^j;,/3^r,/3- (6-4) 

Los terminos de dispersion hacia atras y umklapp son irrelevantes si nos restringimos a 
interacciones repulsivas en el primer caso, y estamos lejos del Uenado medio en el segundo. 

En este capitulo calculamos las funciones de correlacion de los operadores que representan 
ondas de densidad de carga (CDW) y spin (SDW), y superconductividad de tipo singulete 
(SS) y triplete (TS) en el modelo anteriormente presentado. Las funciones de correlacion 
para estos operadores son bien conocidas en ausencia de acoplamiento SO [17,48,50], in- 
cluyendo factores de correcion logaritmica que se originan en terminos irrelevantes [85,115] y 
dependencia con el tiempo y la temperatura [48] . En el presente capi'tulo se extienden estos 
calculos al caso en que las interacciones SO se encuentran presentes, y se estudia como se 
modifican los exponentes de sus decaimientos algebraicos. Como resultado encontramos in- 
teresantes modificaciones del diagrama de fases del sistema. Para ciertas regiones del espacio 
de parametros, la interaccion SO cambia la fuctuacion dominante, y hace que susceptibili- 
dades que eran finitas, ahora se vuelvan tambien divergentes a bajas temperaturas. 

6.2.2 Bosonizacion 

El Hamiltoniano ()(j.2|l puede estudiarse mediante la tecnica de bosonizacion, como en la Ref. 
37,38. Aunque es indistinto emplear bosonizacion funcional u operacional, elegiremos esta 
ultima para Uevar un paralelismo lo mas estrecho posible con dicha referenda, aunque final- 
mente los valores medios seran calculados utilizando tecnicas funcionales. Por conveniencia 
definimos una velocidad promedio vq = {vi + f2)/2 y la diferencia 6v = V2 — fi- A partir 
de la figura Ifj.ll se observa que los momentos de Fermi tambien son diferentes para ambas 
ramas, por esa razon, hacemos lo mismo para los momentos de Fermi, definiendo el promedio 
^0 = (^1 + ^2)72 y la diferencia Sk = k2 — ki. Para efectuar la bosonizacion introducimos 
los campos de fase usuales (pp y 0o- para los grados de libertad de carga y spin, y los corre- 
spondientes campos duales Up y 11^^ siguiendo el procedimiento explicado en el capitulo El 
En terminos de los campos bosonicos, el Hamiltoniano puede representarse como 
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+ 5vjdx [(9,0p)n. + {dM^pl ■ (6.5) 

Vp^a son las velocidades de propagacion de los modeos colectivos del modelo desacoplado 
{6v = 0), y Kp f^ son las constantes de dureza. El acoplamiento SO aparece como un efecto 
que rompe la separacion spin-carga, lo cual se manifiesta en la presencia del tercer termino en 
la ultima ecuacion. Sin embargo, este Hamiltoniano puede ser diagonalizado en terminos de 
dos nuevos campos de fase, que portan una mezcla de carga y spin. Dado que el Hamiltoniano 
posee terminos cruzados en los campos y los momentos, su diagonalizacion no es trivial, por 
ejemplo no puede diagonalizarse por una transformacion de similitud, porque resulta no 
canonica. Dejamos para el apendice iB] los detalles. Baste mencionar que las velocidades de 
propagacion de estos nuevos modos colectivos son 





K K 



p/ J 



(6.6) 



A medida que 5v ^ max(fp,fo-) y "W- — > mmivp^Va). A medida que 5v se 

incrementa, f_ disminuye hasta anularse en los puntos 



Kp 



(6.7) 
(6.8) 



En estos puntos, el congelamiento del modo bosonico mas lento esta acompafiado por una 
divergencia en las funciones respuesta de carga y spin. La compresibilidad de carga estatica 
K diverge para 6v = 6vp y para 6v = 5Vfj ocurre una divergencia en la susceptibilidad de spin 
estatica x- Sus comportamientos son 



K =Kq 



6v 
6vp 
6v 
5v„ 



Xo 



2Kp 

TTVp 

2K„ 



(6.9) 
(6.10) 



donde y Xo son los valores de k y x en ausencia de acoplamiento SO. Mas alia de estos 
puntos las susceptibilidades se vuelven negativas. Este comportamiento de las funciones 
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respuesta estaticas junto con el hecho de que las velocidades de los modos colectivos se 
anulen indican que el sistema se ha vuelto inestable [17, 116] y que realiza una transicion 
de fase de primer orden [85]. Para Kp > K^, Svp resulta menor que 6va, y a medida que 
6v crece desde 0, la divergencia fisica tiene lugar en la compresibilidad de carga. Esta 
inestabilidad se conoce como separacion de fases y se ha mostrado que ocurre en el modelo 
de Hubbard extendido [117] y en el modelo t — J [118]. En el caso en que Kp < K„, la 
inestabilidad tiene lugar en el subsistema de spin, y esta relacionada a la llamada transicion 
metamagnetica, observada por ejemplo en el compuesto cuasi unidimensional Ba3Cu204Cl2 
[119]. Tambien ocurre en el diagrama de fases del modelo XXZ con segundos vecinos [120]. 
En presencia de un potencial quimico (campo magnetico), la region en la cual k (x) se vuelve 
negativa esta asociada a la coexistencia de dos fases con diferente concentracion de agujeros 
(magnetizacion) . La divergencia en k fue hallada tambien en otros modelos con dispersion 
asimetrica y sin interacciones SO [46]. 

6.3 Funciones de correlacion 

Enfoquemos ahora nuestra atencion en las funciones de correlacion. Los operadores de interes 
son 

0CDW 
0SDW,x 
0SDW,j/ 
0SDW,z 

0SS 

0TS,O 
0TS,a 

Para obtener su forma bosonica basta utilizar las expresiones de equivalencia entre cam- 
pos fermionicos y bosonicos, ()2.60|1 y ()2.76|1 . El resultado es 



r,a 
r,a 

r,a 

^ a 



-2irkpx 



Ra- 



(6.11) 
(6.12) 

(6.13) 
(6.14) 

(6.15) 

(6.16) 
(6.17) 
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2 I — I — 

0CDW = — cos(2A;oa; + v27r</)p) COS V 27100-, (6.18) 
7ra 

2 / — 
0SDW,x = — cos(2/coa; + v 27r0p) cos(5A;x + V^nO^j), (6.19) 
na 

2 / — 
0SDW,y = — cos(2A;oa^ + v 27r0p) sin((5A;x + V2n9t^), (6.20) 

2 

0SDW,2 = — sin(2fcoa; + v27r0p) sin v 27100, (6.21) 
7ra 



0SS = — — e-^^^" sin v^0., (6.22) 
V27ra 

0Ts,o = -^e"*^^" cos v^0o, (6.23) 
V27ra 

27ra 

Recordamos que, segiin se definio en el capituloEl Oi, esta vinculado con por la relacion 

6.3.1 Funciones de correlacion a temperatura finita 

Las funciones de correlacion se calcularon en el marco de la integral funcional dentro del 
formalismo de tiempo imaginario de Matsubara [47]. En este formalismo, se definen como 



i?,(x,r;/3) = (0,(a;,r)0l(O,O)^ 

= T^^pT^^a'Dct>pct>^&,{x, r)0l(O, 0) exp{-^[n„ 0,]}, (6.25) 

donde Zq es la funcion de particion 

Zo = y"pnpPnoP0pP0oexp{-S[n„0,]}, (6.26) 



S es la accion euclidea 



rP fP r 

S[U^,(f)^]= dTn{T)-i dr dxIi^{x,T)d^(j)^{x,T), (6.27) 
Jo io J 

y r es el tiempo imaginario. Las propiedades de tiempo real se obtienen por continuacion 

analitica r — > it. Como los operadores ()6.18|) estan expresados en terminos de (pi, y 6^, es 

conveniente trabajar en terminos de Oy en lugar de 11,^. La accion expresada en las nuevas 

variables es 
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/3 



dr I dx 



+ 



dr / dx 



+ 



dr / dx [djp {idr(t)p + 6vd^(j)^) + dj^ {idr(t)a + 5w<9x0p)] • (6.28) 



Por la forma que poseen los operadores (jfi.lHj) . sus funciones de correlacion pueden es- 
cribirse de forma general como una combinacion de terminos de la forma 



exp O ^ 13k [v? 



k[x) - (pk 



(6.29) 



donde definimos el campo if de la siguiente manera: 



(6.30) 



con l3k constantes apropiadas, y j = 1, 2, 3, 4. La expresion ()6.29|) es igual a (formula obtenida 
en el apendice EI) 



expJ^A/?, [A^-'(:^)-A,^.i(0)] i, 



(6.31) 



donde 



(^,(x)^,(y)) = A-.i(x-y). (6.32) 

De modo tal que las funciones de correlacion que nos interesan (exponenciales de los 
campos) quedan expresadas en terminos de funciones de correlacion de los campos (f)u J 
6y. Estas se pueden calcular siguiendo el procedimiento estandar de construir una funcional 
generatriz con la accion ()(j.28|l . y derivar funcionalmente respecto de las fuentes externas. 



Mediante este procedimiento se encuentra 



n^,,^^{x,T-[3) = {[iPi{x) - ^i(O)] [ipj{x) - iPj{Q)]) 

= ;^E / dk{l-e-''^-'--^)G^^,^^{uJn,k)e-^\^\. (6.33) 

^ n=-oo -^^^ 
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donde cUn son las frecuencias de Matsubara cOn = 27m/ p. El factor e '"''^l actiia como regulador 
en el U.V. Las funciones G^^^^.{un, k), simetricas en el intercambio de ipi con ipj, resultan 

KpVp{ujl + k'^vl) - SvK^Vak'^ 



G(f) 

Grp 

Ge 
G^ 

Gj, 



[UJr 



[UJr 



k) 



k) 
k) 
k) 

k) =i 
k) =i 





+ vlk^){ul + vlk^) 




uol + k'^v'l) - dvKpVpk'^ 




+ vlk^){uol + vlk^) 




+ k'^^l) - SvKpV^k'^ 






_ KpVa{ 


ujI + k'^vl) - SvK^Vpk"^ 



k) =6v 
k) =5v 
k) =6v- 
k) =6v 



KpKAujl + vlk-'){ujl + vlk^) 

UOn [uol + k^ {vl + ^V'^)] 

k{ujl + vlk^){ujl + vlk^) 

[ul + e {vl + 5v^)] 
k{ujl + vlk^){ujl + vlk^) 
iuJnk{KpVp + K„v„) 



{ujl + vlk-^){ujl + vlk^) 
iuJrik{KpV^ + K^Vp) 
KpKAujl + vlk^){ul + vlk^) 
e{KpVpV^-KJv^)/K^-ujl 



{ul + vlk^){ujl + vH^) 



k\K„v,Vp - Kp5v^)/Kp 



{ujl + vlk^){ujl + vlk^) ■ 

Debe notarse que las funciones que involucran mezclas de campos de carga con campos de 
spin son proporcionales a 6v, de modo que se anulan cuando no hay acoplamiento spin-orbita, 
restaurando la separacion spin-carga. Dejaremos los detalles del calculo de las integrales para 
el apendice|Xl y aqui presentaremos los resultados finales para las funciones ()(i.25j) : 



^CDw(a;, T] (3) = RsBw,z{^^ /5) 



cos2/cox 
2(7ra)2 



2:4-2:4 



-^(Kpp'L+K„u'^_)/2 



H sign(a;T) 



+ h.c. 



(6.44) 



2(7ra)2 
cos 2k2X 



+ 



2iTia) 



-(^KpV<',+ii\IKa)l2^e'\ 



{z.-z.) 



-{Kpu''_+^'L/K^)/2+eZ 



(6.45) 
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±ii5fea; /- \ ±Gsign(a;T) 

(zTraj \Z-\-Z- J 

(6.47) 

donde 



sin -^(e + v+T + ix) 

.± = ,f ^ (6.48) 

sm 

sin ^^(e + v+T — ix) 

z+ ^ ^-'\ ,f ^ (6.49) 

sm ^ 

y los exponentes dependen de las constantes K multiplicadas por factores que incluyen 
dependencias en las velocidades. Estan dados por 

vi^± —— ^ ^ — ^ 6.50 

v± vj^ — vt 

ni = ±—— ^ „ ' ^' 6.51 

v± vj^ — vt 

5v vl - (5vl - 5v^) 
9i = ± ± — ^ 6.52 

V± V± — VI 



con X — p,a,y 

H ^5.^^4±^ (6.53) 
v^ —VI 

G =^. ^'-/-^+'y^ . (6.54) 
v^ — VI 

f^. y /i± son positivos, y 6'^, G j H tienen el mismo signo que 6v. En el modelo sin 
acoplamiento SO, la simetria SU(2) puede restaurarse fijando el valor de K^^ — 1, valor que 
emerge naturalmente si el modelo bajo estudio es el Kmite contmuo de un modelo en la red 
con solamente interacciones de tipo densidad de carga. En nuestro caso esta simetria esta 
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explicitamente rota desde el principio, y no es posible restaurarla. Esta ruptura se manifiesta 
en las diferencias en los decaimientos entre las funciones de correlacion de operadores SDW 
en la direccion z y las direcciones x e y. 

Igual que en el caso en que no hay acoplamiento SO, las funciones de correlacion para 
operadores SDW en la direccion z y operadores CDW son iguales. Lo mismo ocurre con 
las funciones de correlacion para operadores TS y SS. Esta degeneracion se rompe al incluir 
correcciones logan'tmicas que surgen si se tienen en cuenta terminos irrelevantes de umklapp 
y backscattering [115]. 

Un punto interesante de observar es la aparicion de dos terminos en la funcion SDW,a;?/ 
[Ec. ()f).45j) ] donde las modulaciones poseen diferentes frecuencias y decaen con diferentes 
exponentes. Como 9^ tiene el mismo signo que 5v [ver la Ec. ()6.52p y el comentario debajo 
de la Ec. fl6.54|) ] para V2 > fi (f2 < "^1) el termino dominante es el de frecuencia ^2 (fci). En 
otras palabras, la mayor frecuencia domina. Ademas -Rts,±i se vuelve oscilante. 

6.3.2 Funciones de correlacion instantaneas a temperatura cero 

Hasta aqui liemos obtenido formulas muy generales para las funciones de correlacion de- 
pendientes del espacio, tiempo imaginario y la temperatura. Podemos lograr una mayor 
comprension de la fisica del problema observando el decaimiento algebraico de las funciones 
de correlacion instantaneas (r = 0) y a temperatura cero y estudiar como sus exponentes se 
modifican respecto del caso con acoplamiento SO nulo. El comportamiento general de estas 
funciones es 

Ri{x) ~ \x\-^+"\ (6.55) 

Los exponentes a^s determinan la divergencia de la correspondiente susceptibilidad en el 
espacio de Fourier cuando T — > 0, Xi(T) ~ T~"' [17]. De este modo, estas inestabilidades 
resultan de una naturaleza completamente diferente a las descriptas en las Ecs. ()6.9|) y 
fl6.10|) . Las expresiones obtenidas para los ctj son 



acDW = asDW,2 =2 - KpU^ - K^u"" (6.56) 

«SDW,. = «SDW,y =2(1 + |^"|) - K.iyP - fi'^/K, (6.57) 

«ss = axs.o =2(1 + \ef\) - ^iP/Kp - K^v" (6.58) 

aTS,±i =2 - ^l'/Kp - fi^K^. (6.59) 



Estos son los nuevos exponentes, que retienen la misma estructura que en el caso de SO nulo, 
pero modificados por los factores 
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Figura 6.2: Comportamiento de los exponentes a^s como funcion de 5v (en unidades de 
Vq). Para Vp = 1.2vo, = O.Svq, Kp = 0.6 y K^^ = 0.85. Para 6v > 0.16 las fluctuaciones 
SDW,xy se vuelven dominantes, y para 6v > 0.25 ass se vuelve positive, y xss divergente 
para T — ^ 0. 



/ = /i^ + /i^ (6.60) 
z/^ = z/^ + (6.61) 
e^ = 9l + 9^. (6.62) 

Cuando 5f — 0, se verifica que 6''^ 0, y /i'^, i/'^ 1, de modo que reproducimos los 
resultados correctos para el caso SO nulo. 

Para acoplamiento SO finito, 6v es un parametro que juega un rol en determinar cual 
es la funcion de correlacion que decae mas lentamente y cuales son las susceptibilidades 
divergentes. En la Fig. 16.21 observamos, a modo de ejemplo, el comportamiento de los 
exponentes como funcion de 6v para Vp = 1.2vq, = O.Svq, Kp = 0.6 y = 0.85. Para 
Sv pequeiio, las fluctuaciones CDW son dominantes, pero para 6v > 0.16^0 las correlaciones 
SDW,xy decaen mas lentamente. Para 6v pequeno, las fluctuaciones CDW y SDW son las 
linicas fluctuaciones divergentes para T — > 0, pero para 6v > 0.25fo, ass se vuelve positivo 
y xss divergente para T — >• 0. Calculos de la estructura de bandas electronicas modificadas 



76 



Kp 




Figura 6.3: Diagrama de fases en el espacio Kp — K„. Las fases entre corchetes son 
dominantes, que se vuelven dominantes para acoplamiento SO suficientemente fuerte. 
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Figura 6.4: Diagrama de fases en el espacio Kp — 5v para Vp = 1.2vq, v„ = O.Svq y diferentes 
valores de K^. 6v > 5v„ debajo de la Imea punteada y ocurre el metamagnetismo. 
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por el acoplamiento SO muestran que estos valores de 6v deberian corresponder a sistemas 
de electrones cuasi unidimensionales ti'picos. 

Un analisis cuidadoso de los exponentes nos permite construir un diagrama de fases en 
el espacio Kp — K„. (Fig. lfi.H|l . En cada region indicamos las fluctuaciones dominantes 
para 6v pequeiio, y entre corchetes las dominantes para Sv mayores. Otras fluctuaciones 
sub dominantes no estan indicadas. En la Fig. 16.41 se observan cortes del diagrama de 
fases en el espacio Kp — 6v para Kp < 1 y diferentes valores de K^^. Para 6v pequeiio las 
fluctuaciones CDW dominan, y para 6v mas grande , el sistema puede encontrarse en la 
fase SDW o SS dependiendo de los valores de Kp y K^. En la region debajo de la Ifnea 
punteada, 6v < 6va, la susceptibilidad estatica de spin se vuelve negativa, y tiene lugar el 
metamagnetismo. 

6.4 Conclusiones 

En este capftulo hemos calculado funciones de correlacion para operadores de fluctuaciones 
de ondas de densidad de carga y spin, y superconductividad singulete y triplete en un mod- 
elo de electrones altamente correlacionados en una dimension con acoplamiento spin-orbita. 
El calculo se realizo en funcion de la temperatura, y al final se estudiaron las funciones 
instantaneas y a temperatura cero. El acoplamiento spin-orbita destruye la separacion spin- 
carga como se mostro en las Ref. 37, 38 y modifica los exponentes de los decaimientos de 
las correlaciones. Como consecuencia se modifica el diagrama de fases del sistema. Un 
acoplamiento spin-orbita suficientemente fuerte, es responsable de un cambio en las fluctua- 
ciones dominantes y de promover que nuevas susceptibilidades se vuelvan divergentes para 
T ^ 0. 
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Capitulo 7 
Conclusiones 



En esta tesis hemos estudiado diversos aspectos teoricos de los sistemas de electrones al- 
tamente correlacionados en una dimension espacial. En particular hicimos hincapie en el 
tratamiento del modelo de Tomonaga-Luttinger con spin y algunas de sus extensiones, me- 
diante el procedimiento de la bosonizacion. Al aplicar este metodo a sistemas de materia 
condensada, el primer punto en el que hay que detenerse es en el calculo del determinante 
fermionico asociado a la teorfa. Como el determinante esta mal definido se debe implemen- 
tar un mecanismo de regularizacion. El metodo elegido viene usualmente dictado por las 
simetrias que posea el modelo en cuestion. En teon'as de materia condensada la ausencia 
de invarianza de Lorentz, otorga una libertad aiin mayor para elegir el regulador. En el 
capitulo El mostramos que al utilizar el regulador usual que preserva la invarianza de Gauge 
y de Lorentz en teorias de campos relativistas se arriba a resultados que no concuerdan con 
los obtenidos mediante el camino usual de bosonizacion operacional en materia condensada, 
y encontramos la forma exacta que debe tener este regulador para reproducirlos. Como coro- 
lario obtenemos el jacobiano de las transformaciones quirales, base del desacople fermionico 
mediante integrales funcionales [41]. Desafortunadamente, un principio fisico que sirva de 
gui'a para elegir a priori el regulador todavi'a falta, aunque lo mismo ocurre en el enfoque 
operacional. 

En el capitulo El empleamos el metodo de bosonizacion funcional, mejorado de acuerdo 
a lo diclio anteriormente, para estudiar el modelo de Tliirring no local con interacciones 
que invierten el spin electronico (los tratamientos previos de las interacciones de inversion 
de spin condujeron a un modelo no abeliano en el cual el analisis del contenido fisico se 
vuelve engorroso). Restringiendonos al caso de interacciones de inversion de spin locales, 
obtuvimos una accion bosonica efectiva cuyos grados de libertad de carga coinciden con los 
encontrados en el capitulo El al tratar el problema de dispersion hacia adelante sin spin. 
Respecto del sector de spin, que es el de mayor interes, liallamos que se corresponde con un 
modelo seno-Gordon no local, cuya integrabilidad, en contraste con el caso local, no ha sido 
demostrada hasta el momento. A pesar de este hecho, pudimos mostrar explfcitamente que 
ambos sectores se desacoplan dando lugar a la separacion spin-carga [42]. 
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En el capitulo El revemos la aproximacion armonica autoconsistente, y su aplicacion a 
sistemas de materia condensada. En particular hacemos una extension que permite atacar 
problemas no locales, y obtenemos una expresion para el gap del sector de spin del modelo 
introducido en el capitulo lUconio funcion de los potenciales de dispersion hacia adelante. Por 
otro lado, como una leve digresion con respecto a la linea principal de esta tesis, proponemos 
un camino diferente al usual para la determinacion del parametro asociado a la aproximacion, 
basado en una consecuencia del teorema c de Zamolodchikov en teon'as de campos conformes. 
Verificamos la validez de la nueva tecnica en el modelo seno-Gordon, en el que obtenemos 
una mejora en los valores aproximados para la masa del soliton como funcion de la constante 
de acoplamiento /5 del modelo, con respecto a los predichos mediante el procedimiento usual. 
Por ultimo aplicamos la nueva tecnica al estudio de un problema no trivial, para el cual 
existen escasos resultados analiticos: el modelo de Ising bidimensional a, T ^ T^. y h ^ {). 
Hallamos una expresion aproximada que relaciona la longitud de correlacion C,, T — Tcj h j 
que resolvemos numericamente para obtener ^ como funcion de T — Tc j h. [43] . 

Finalmente en el capitulo El calculamos funciones de correlacion en sistemas de electrones 
altamente correlacionados en una dimension en los que los grados de libertad de carga y spin 
se encuentran acoplados a traves de la interaccion spin-orbita. Esta interaccion rompe la sep- 
aracion spin-carga [37, 38] y modifica los exponentes de los decaimientos de las correlaciones. 
Estudiamos fluctuaciones de tipo ondas de densidad de carga y spin, y de tipo supercon- 
ductor singulete y triplete como funcion del espacio-tiempo euclideo y de la temperatura. 
Ademas investigamos las funciones de correlacion instantaneas a temperatura cero, que nos 
permiten extraer los exponentes criticos. Mostramos que la interaccion spin-orbita modifica 
los exponentes del decaimiento de las funciones de correlacion y el diagrama de fases del 
sistema. Ademas encontramos que susceptibilidades que eran finitas a bajas temperaturas, 
pueden aliora volverse divergentes de la interaccion spin-orbita [44]. 
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Apendice A 

Propagador bosonico a temperatura 
finita 

En este apendice describiremos el calculo del propagador [16,39,48] 

A \x,r;(3) = - / --^--^ (A.l) 



(3 ^ J_^27ru;l + v^k^ 

' n=—oo °" " 

donde Un — 2n7r//3 son las frecuencias de Matsubara y r e [0, /3] (A~^(x, t;/3) es periodica 
en r con pen'odo P). Esta funcion es divergente en el infrarojo, pero nosotros en realidad 
estamos interesados en la combinacion 

-1 oo „oo J7. -ikx-iu!„T _ 1 

n=— oo " 'J" " 

que es regular en esa region. El termino sustraido, sin embargo, es divergente en el ultravi- 
oleta, razon por la cual hemos agregado un factor e"^''^' que regulariza dicha divergencia. 

A continuacion escribimos e""^""^ en terminos del seno y el coseno y nos quedamos 
linicamente con el termino del coseno, ya que el del seno es nulo por simetrfa. La sumatoria 
resultante puede ser facilmente evaluada mediante el metodo de los residuos, o consultada 
laRef. 121: 



cosna; 7rcoslia(7r — a;) 1 /a q\ 

^ -n? + a? 2a sinh an 2a^ 

n=l 

Utilizando este resultado, obtenemos entonces 

^ ^ ^ ' J-oo 47rt;|/c| sinh^ ^ ' 
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Para calcular la integral antes debemos partir la integral en y en el termino integrado 
en la region negativa de k hacer el cambio k — ^ —k. Ademas reescribimos las funciones 
hiperbolicas en terminos de exponenciales, obteniendo 







r°° rlh (p-ikx _|_ ikx\ ( p-kvT i -kv((3-T)\ _ O /' 1 i p~kv(3\ 

(A.5) 

e-^ ^' \ ^ + (s - s*) (A.6) 



Anvk 1 — e ^'^ 

oo 



Anvk 1 — 



(A.7) 



donde s = f r + ia: y r = u/?. Esta ultima integral tiene forma estandar, (3.413(.l)) en la 
Ref. 121. Con este resultado, 

donde F es la funcion gamma. Para simplificar este resultado reemplazamos en la ecuacion 
anterior 1 + e/r por 1 — e/r, cambio chico si e es chico, ademas usamos que r(z)r(l — z) = 
7t/ siuTTZ para dar finalmente. 

A-\x, T- (3) - A-i(0, 0; f3) ~ In . '^''^l'' / ^ + ^ (^'9) 

47rf sm 7r(e/r + s/r) 

X sin ■^'^ 



^ In—— ^^ + (c.c.). (A.IO) 

Ativ sin ^{e + vt + ix) 

Las restantes integrales ()(i.33|l - ()fi.43|l se evaluan utilizando el metodo de fracciones sim- 
ples, seguido del procedimiento descripto en este apendice. A continuacion damos los resul- 
tados: 



oo 



ii=^ r dk ^ 



(1 - e"**^^-^'^"") e 



E|fc|. 



n=— oo 
^2" 



27r6 \ \?-c^ 



In 



sin -^{e + br 



ujl + b^k^){iul + c^k^) 

1 2 



^2;) 1 f — sin 7^ (e + cr + «x) 



+ 



sm 



6/3 



27rc 
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In- 



sm 



c/3 



+ C.C. 

(A.ll) 
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-| /"OO 



-ikx—iu}nT\ 



-e\k\_ 



Sign XT 
27r 



n=— OO 
1 
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,.7r l .TT 

m sm 77;^ (e + or — za; j + — — m sm — (e + cr — za; j 



6/3 



c/3^ 



+ C.C. (A.12) 



1 /"^^ °° 



sign XT 
27r 



-ifca;— ia/nTj 
TT 



6 — d 

— ) In sin --(e + 6t - ix) + 

— J bp 



k{u;l + b^k^){ujl + c^k^) 

n 



c^-b^ 



In sin —(e + cr — ix) 



+ c.c. 
(A.13) 
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Apendice B 



Diagonalizacion del Hamiltoniano 
bosonico 



Supongamos por simplicidad, una version de dos grados de libertad de nuestro Hamiltoniano 
bosonico 16.51 de la forma 



Para diagonalizar este Hamiltoniano no podemos emplear una transformacion de similitud 
porque resulta no canonica debido a los terminos que mezclan coordenadas e impulsos. En 
su lugar emplearemos metodos simplecticos como se describen en la segunda edicion del libro 
de Goldstein [122]. 

En forma general, tomamos un Hamiltoniano de n grados de libertad y construimos un 
vector ^ con 2n elementos compuesto por las n coordenadas qi y los n momentos pj. Las 
ecuaciones de Hamilton en notacion simplectica se escriben 



donde J, que es la matriz cuadrada de 2n x 2n compuesta por las matrices de n x n nula e 
identidad, segun el esquema 



Mencionamos algunas propiedades de estas matrices. Su cuadrado es la identidad cam- 
biada de signo 






(B.3) 



P = -1. 



(B.4) 
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Es tambien ortogonal: 

JJ = 1 (B.5) 

con lo cual 

J = -J = r\ (B.6) 

De su ortogonalidad se deduce que el cuadrado de su determinante es 1, si bien se puede 
probar la afirmacion mas fuerte 

detJ = +l. (B.7) 

Una transformacion canonica es una tranformacion de las coordenadas y los impulsos de 
la forma 



Qi^Qi{q,p), Pi^Pi{p,q) (B.8) 

donde Pi y Qi son las nuevas coordenadas e impulsos, que preserva sus conmutadores 
cuanticos (o corchetes de Poisson en cl caso clasico). En notacion simplectica, si defini- 
mos al vector r) de las nuevas coordenadas e impulsos, la transformacion resulta 

e = e(^), (B.9) 

y la condicion para que la transformacion resulte canonica (Uamada condicion simplectica) 
se escribe 

MJM = J (B.IO) 

donde M es la matriz jacobiana de la transformacion y M su matriz traspuesta. 
Sea H un Hamiltoniano cuadratico general de n grados de libertad 

H^l^SC (B.ll) 

donde S es una matriz cuadrada simetrica constante. Para hallar la transformacion canonica 
lineal que lo diagonaliza, apelaremos a las ecuaciones clasicas de movimiento. Las ecuaciones 
de Hamilton para este sistema resultan 

e = JSe. (B.12) 
Para resolverlas se puedc proponer el ansatz 

C = ae^'"* (B.13) 
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donde a es un vector constante. La ecuacion para a resultante de reemplazar este ansatz en 
f|RT2|l es 



-iJSa = ua (B.14) 

es una ecuacion de autovectores, donde los autovalores u resultan las frecuencias carac- 
teristicas del sistema. Sea U una matriz de un conjunto posible de vectores propios de 
— zJS ordenados en columnas tales que solo sus direcciones estan fijas, no sus modulos, y sea 
D = diag(ci;i, ...,ujn, —^i, —^^n) la niatriz de los autovalores. La ecuacion ()B.14|) se escribe 
en terminos de U y D como 



-zJSU = UD. (B.15) 

Ahora bien, la matriz U no es en general la matriz de ninguna transformacion canonica, ya 
que no cumple la condicion simplectica ()B.10|) . Mostraremos que es posible fijar condiciones 
sobre los modulos Ck de los autovectores de modo tal que la nueva matriz de autovectores 
cumpla con dicha condicion. Para ello introducimos la matriz de los modulos C = diagc^, 
y definimos M = CU como una nueva matriz de autovectores, donde cada autovector esta 
multiplicado por una constante a determinar. Si le exijimos a M que cumpla la condicion 
simplectica, entonces 

OCJCU = J, (B.16) 
y multiplicando por U^^ a izquierda y por U"^ a derecha obtenemos 

CJC = U-^JU-^ (B.17) 

que resultan ecuaciones cuadraticas para los c^. 

Una vez obtenida la matriz de la transformacion canonica, podemos transformar el Hamil- 
toniano reemplazando ^ = M?7, lo que resulta 

H = l^S^ = lvMSM7^. (B.18) 
Por la propiedad ()B.4|1 . H se puede escribir 



=_!r/MJ(-zJS)M?7 (B.19) 
= ^^MJMD77 (B.20) 
= Y^JD^, (B.21) 
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donde hemos usado la ecuacion de autovectores ()B.15|) . que tambien satisface M, y la condi- 
cion simplectica para M. Finalmente si empleamos la definicion de D y J, en terminos de las 
nuevas coordenadas e impulsos qi y pi (que componen rj), el Hamiltoniano se escribe 

H = iujiqipi (B.22) 
Por ultimo, mediante la transformacion canonica 



qi=Q^-^ — , (B.23) 



UJi 



2 2 

H resulta diagonal: 



Pi =17 - (B.24) 



H = P! + ^Ql (B.25) 



2 
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Apendice C 

Identidad litil para el calculo de 
valores medios bosonicos 



En este apendice probaremos la identidad 



exp IzJ^Pk [Mx) - MO)] ) = exp i J] A/?, [Ari(x) - A-.^(O)] \ [C.l] 



donde (/? es un campo bosonico y A.^.^(x) es el valor medio bosonico 



La identidad es valida para valores medios calculados con acciones bosonicas cuadraticas 
en el campo ip. Sea S una accion de este tipo 



^ =\ j d^xipi{x)Aijipj{x) (C.3) 

donde Ajj es un operador simetrico en los indices i y j- La funcion A~.^(x) resulta ser 
entonces la funcion de Green de este operador. El lado izquierdo de ()C.1|) se puede calcular 
mediante el siguiente procedimiento: 
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exp <^ I 

k 



Y,Pk[M^)-Mo)]\)= (C.4) 



= "Cv^expl-i J d'^x' (pi{x')Aij (pj{x') + i^/3k[(pkix) - (pkiO)]j (C.5) 

= ^ J ^V^expj-i J (fx'ipi{x')Aijipj{x') + j (fx' ipk{x')i[3k[5{x' - x) - 5{x')] 

(C.6) 

= j "^V^expj-i j (fx'ipi{x')Aij(pj{x')-j (fx'ipk{x')jk{x',x)^ (C.7) 
donde Zq es la funcion de particion 



Zo = J P<^exp|-^ J d^x'ipi{x')Aijip{x')^ (C.8) 



y definimos 



jk{x',x) = -if3k [S{x - x) - 8{x')] . (C.9) 
que es independiente de los campos. A continuacion hacemos una traslacion en el campo ip 

^i(x') ^ + j d^yA-\x' - y)jkiy,x). (C.IO) 

Este es el unico cambio en la integral funcional porque la medida de integracion funcional 
es invariante frente a traslaciones. Finalmente, el valor medio queda 



exp yY^f^k ['^'^(^) - (C.ll) 

= exp|^ J d^x' d%j,{x',x)Ar/{x' -y)j,{y,x)^ (C.12) 

= [^^'^^ ^^^^^"^^ " } ^^-^^^ 

donde en el ultimo paso hemos reemplazado jk{x',x) por su definicion ()C.9|) . e integramos 
en x' e y. 
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